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In this review, we summarize tests of standard electroweak (EW) theory at the highest available 
energies as a precursor to the Large Hadron Collider (LHC) era. Our primary focus is on the 
published results (as of March 2010) from proton-antiproton collisions at -y^ = 1.96 TeV at the 
Fermilab Tevatron collected using the CDF and DO detectors. This review is very timely since 
the LHC scientific program is nearly underway with the first high-energy (y/s = 7 TeV) collisions 
about to begin. After presenting an overview of the EW sector of the standard model, we provide 
a summary of current experimental tests of EW theory. These include gauge boson properties 
and self-couplings, tests of EW physics from top quark sector, and searches for the Higgs boson. 
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I. INTRODUCTION 

The goal of particle physics is to explain the nature of 
the Universe at its most fundamental level, including the 
basic constituents of matter and their interactions. The 
standard model of particle physics (SM) is a quantum 
field theory based on the SU{?>)c'^SU{2)i,®U[1)y gauge 
symmetry group which describes the strong, weak, and 
electromagnetic interactions among fundamental parti- 
cles. This theory has been the focus of intense scrutiny 
by experimental physicists, most notably at high energy 
particle colliders, over the last three decades. It has been 
demonstrated to accurately describe fundamental parti- 
cles and their interactions up to O(IOO) GeV, with the 
existence of non-zero neutrino masses and mixing being 
the only known exception. Despite the success of the SM, 
there are many reasons to believe that the SM is an effec- 
tive theory which is only valid up to ~1 TeV. Some of the 
shortcomings of the SM will be described in Section [TTI 

Particle physics is embarking on a unique, and possibly 
defining, period in its history with the start of particle 
collisions at the Large Hadron Collider (LHC) at CERN. 
At the LHC, bunches of protons will be collided with a 
planned 14 TeV of center-of-mass energy, creating condi- 
tions that existed only a tiny fraction of a second after 
the big bang. This is seven times the center of mass 
energy of collisions at the Fermilab Tevatron. For the 
first time, physicists will be able to directly probe the 
TeV energy scale in the laboratory, where new physics 
beyond the SM, with the potential to revolutionize our 
understanding of the Universe, could be apparent. We 
can only speculate about what form this will take. Is is a 
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tantalizing prospect that on the horizon is a revolution in 
our understanding of the Universe that includes a more 
complete theory of particles and their interactions which 
may explain dark matter and dark energy. 

The purpose of this review is to present tests of the 
electroweak (EW) sector of the SM (S'[/(2)l «> C/(l)y) at 
the highest available energies as a precursor to the LHC 
era. Our focus is on published results from collider data 
collected using the CDF ^ and DO [l] detectors during 
Run II at the Fermilab Tevatron as it relates to our un- 
derstanding of electroweak interactions and spontaneous 
symmetry breaking. After an overview of electroweak 
theory (Section [ll]), we present current results on gauge 
boson properties and self-couplings (Section Hill) , tests of 
electroweak physics from top quark physics fSectionllVp. 
and searches for the Higgs boson (Section N} . 

II. OVERVIEW 

The SM is an extremely successful theory of the strong, 
weak, and electromagnetic interactions. It is based on 
three generations of quarks and leptons, interacting via 
an SU{3)c <S) 5'J7(2)l (E) U{1)y gauge symmetry The 
S'C/(2)l U{1)y symmetry is spontaneously broken to 
electromagnetism, C/(1)_em, by the vacuum-expectation 
value of the Higgs field. Given this field content and 
gauge symmetry, the most general theory that follows 
from writing down every term of dimension four or less 
is the SM. 

In the SM, as usually understood, neutrinos are ex- 
actly massless, and do not mix. Since neutrino mixing 
has been definitively observed, we must go beyond the 
SM in order to describe this phenomenon. There are two 
ways to do this. One way is to extend the field content 
of the model by adding additional fermion and/or Higgs 
fields (e.g. a right-handed neutrino or a Higgs triplet). 
The other way is to extend the SM by adding operators 
of dimensionality greater than four. There is only one op- 
erator of dimension five allowed by the gauge symmetries 
i, 

£5 = ^{I^'^edp)C{feU) + h.c. (1) 

where is the lepton doublet field of the i*'' genera- 
tion and (j) is the Higgs doublet field (the 2x2 matrix 
e and the 4x4 charge-conjugation matrix C are present 
to ensure invariance under S'C/(2)l and Lorentz transfor- 
mations, respectively). When the Higgs doublet acquires 
a vacuum-expectation value, (0) = (0, w/-\/2) {v — 246 
GeV), this term gives rise to a (Majorana) mass for neu- 
trinos, 

C5 = -^M'Ji^'^Ciy' +h.c. (2) 

where Ml^ — d^v'^/A is the neutrino mass matrix. Due 
to the tiny inferred masses of neutrinos, the scale A lies 



around 10^^ GeV, assuming c'^ is not much less than 
order unity. 

There are other indications that the SM is not a com- 
plete description of nature, most of them related to grav- 
itation and cosmology. Even with massive neutrinos 
included, the SM particles only constitute 4.6% of the 
present universe, with the remainder in mysterious dark 
matter (23%) and dark energy (72%). Neither dark mat- 
ter and nor dark energy are accommodated in the SM. 
There is no adequate mechanism for baryogenesis (the 
observed excess of baryons over antibaryons) or inflation, 
which is the simplest explanation of the observed temper- 
ature fluctuations in the cosmic microwave background. 
The SM also provides no explanation of the strong CP 
problem: the lack of observed CP violation in the strong 
interaction, which is allowed by the SM. 

If physics beyond the SM lies at an energy scale less 
than 1 TeV, then we should be able to observe it directly 
at high-energy colliders. If it lies at a scale greater than 
1 TeV, then we can parametrize its effects via higher- 
dimension operators, suppressed by inverse powers of the 
scale of new physics. A, exactly as in the case of neu- 
trino masses described above. Other than the dimension- 
five operator responsible for neutrino masses, the lowest- 
dimension operators are of dimension six, and are there- 
fore suppressed by two inverse powers of A. If A is of 
order 10^^ GeV, as suggested by neutrino masses, then 
these operators are so suppressed that they are unobserv- 
able, with the possible exception of baryon-number vio- 
lating operators that mediate nucleon decay. However, 
there could be more than one scale of new physics, and 
if this scale is not much greater than 1 TeV, its effects 
could be observable via dimension six operators. Opera- 
tors of dimension greater than six are suppressed by even 
more inverse powers of A and can be neglected. 

There are many dimension-six operators allowed by the 
SM gauge symmetry There are three ways to detect 
the presence of these operators. The first is to observe 
phenomena that are absolutely forbidden (or extremely 
suppressed) in the SM, such as nucleon decay. The sec- 
ond is to make measurements with such great precision 
that the small effects of the dimension-six operators man- 
ifest themselves. The third is to do experiments at such 
high energy, E, that the effects of these operators, of 
order {E/A)'^, become large. If i? > A then one must 
abandon this formalism, because operators of arbitrarily 
high dimensionality become significant; however, the new 
physics should then be directly observable. If no effects 
beyond the SM are observed, then one can place bounds 
on the coefficients of the dimension-six operators, c/A^, 
where c is a dimensionless number. These bounds apply 
only to the product c/A^, not to c and A^ separately; 
in fact, there could even be two different scales of new 
physics involved (A1A2 in place of A^). 

This approach to physics beyond the SM, dubbed an 
effective-field-theory approach has the advantage of 
being model independent. Whatever new physics lies at 
the scale A, it will induce dimension-six operators, whose 
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only dependence on the new physics hes in their coeffi- 
cients, c/A^. Another advantage of this approach is that 
it is universal; it can be applied both to tree-level and 
loop-level processes, and any ultraviolet divergences that 
appear in loop processes can be absorbed into the coeffi- 
cients of the operators. Thus one need not make any ad 
hoc assumptions about how the ultraviolet divergences 
are cut off. This effective-field-theory approach thus pro- 
vides an excellent framework to parametrize physics be- 
yond the SM [5i;li]. 

Hadron colliders contribute to the study of the elec- 
troweak interactions in three distinct ways. Firstly, be- 
cause they operate at the energy frontier, hadron col- 
liders are uniquely suited to searching for the effects of 
dimension six operators that are suppressed by a factor of 
(E/A)^. Secondly, they are able to contribute to the pre- 
cision measurement of a variety of electroweak processes, 
most notably to the measurement of the W boson mass 
and the top quark mass. Thirdly, they are able to search 
for new particles associated with the electroweak interac- 
tions, in particular the Higgs boson. These three virtues 
of hadron colliders will manifest themselves throughout 
this review. 



A. Electroweak Interactions 

The electroweak theory is a spontaneously broken 
gauge theory based on the gauge group SU {2)i, (g) U{1)y- 
There are three parameters that describe the theory: the 
gauge couplings g and g' , and the order parameter of 
spontaneous symmetry breaking, v. In the SM, this order 
parameter is the vacuum expectation value of a funda- 
mental Higgs field. These parameters are not measured 
directly, but rather inferred from precision electroweak 
measurements. The three measurements that are used 
to fix these parameters are the Fermi constant Gp deter- 
mined from the muon lifetime formula; the fine structure 
constant a, determined from a variety of low-energy ex- 
periments; and the Z boson mass Mz- With these three 
inputs, the predictions of all other electroweak processes 
can be calculated, at least at tree level. 

The level of precision of electroweak measurements is 
such that a tree-level analysis is insufficient, and one must 
go to at least one loop. At this level, one finds that 
predictions depend also on the top quark mass and the 
Higgs boson mass, since these particles appear in loops. 
In fact, a range for the top quark mass was correctly 
predicted by precision electroweak data before the top 
quark was discovered, and the measured mass falls into 
this range. We are now following the same tack with 
the Higgs boson. Remarkably, the precision electroweak 
data imply that the Higgs boson mass is not far above the 
experimental lower bound of ttih > 114 GeV/c^, which 
means that it may be accessible at the Tevatron as well 
as the LHC. 

The electroweak interaction has many other param- 
eters as well. Along with the top quark mass, there 



are the masses of all the other quarks and leptons, as 
well as the elements of the Cabbibo-Kobayashi-Maskawa 
(CKM) quark-mixing matrix and the Maki-Nakagawa- 
Sakata (MNS) lepton-mixing matrix. Most of these mix- 
ing parameters are not measured at the energy frontier, 
with one exception: the CKM element Vtb that describes 
the coupling oi a W boson to a top and bottom quark. 
The only direct measurement of this parameter comes 
from electroweak production of the top quark at hadron 
colliders via a process known as single-top production, 
discussed in Section HV. CI 



B. Electroweak Symmetry Breaking 

The strong and electroweak forces are gauge theories, 
based on the groups SU{3)c and U{1)em, respectively. 
The associated gauge bosons, the gluon and the photon, 
are massless as a consequence of the gauge symmetry. 
We know that the interactions of electroweak bosons with 
fermions as well as with themselves are also governed by 
a gauge theory, with gauge group SU{2)i,(^U{1)y- Why, 
then, are the electroweak bosons, and Z, not mass- 
less, as would be expected of gauge bosons? In the SM, 
the answer is that the electroweak symmetry is sponta- 
neously broken, and that the electroweak gauge bosons 
acquire mass through the Higgs mechanism. This is the 
most plausible explanation of why the interactions ap- 
pear to be those of a gauge theory, despite the fact that 
the gauge bosons are not massless. However, this argu- 
ment leaves completely open the question of how (and 
why) the electroweak symmetry is broken. 

The simplest model of electroweak symmetry breaking 
(EWSB), which is also the original proposal, is based on 
a fundamental scalar field that is an electroweak dou- 
blet carrying hypercharge Y = 1/2. The potential for 
this scalar field is chosen such that its minimum is at 
nonzero field strength. This breaks the electroweak sym- 
metry to U{1)em, as desired. This simple model, which 
can be criticized on several grounds, has withstood the 
test of time. It predicts that there is a scalar particle, 
dubbed the Higgs boson, of unknown mass but with def- 
inite couplings to other particles. The discovery of this 
Higgs particle is one of the driving ambitions of particle 
physicists, and was a primary motivation for the LHC. 

As mentioned in the previous section, this simple 
model is consistent with precision electroweak data with 
a Higgs particle close to the present experimental lower 
bound of 771 if > 114 GeV/c^. This consistency does 
not rule out more exotic possibilities, however, such as 
two (or more) Higgs doublets, Higgs singlets and triplets, 
composite Higgs bosons, and other alternative models of 
electroweak symmetry breaking. 
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III. ELECTROWEAK GAUGE BOSONS 

In the SM, the W and Z bosons mediate the weak 
force and acquire mass through the Higgs mechanism, 
as described in Section |lTl The W boson was discov- 
ered in 1983 in pp colhsions at the CERN SPS by the 
UAl and UA2 experiments B; Hi , with discovery of the 
Z boson soon to foUow [§; The discovery of these 
gauge bosons at CERN represents a dramatic vaUdation 
of Glashow-Salam- Weinberg (GSW) theory which pre- 
dicted the existence of neutral currents mediated by a 
new gauge boson, the Z boson, and predicted the W 
bosons to describe nuclear /3-decay and. Together with 
the massless photon, these comprise the gauge bosons 
of the electroweak interaction. High precision studies of 
the Z boson properties made by the LEP collaborations 
and the SLD collaboration 11] using e+e~ collisions have 
provided stringent tests of electroweak theory. 

The W and Z bosons are copiously produced in pp 
collisions at the Fermilab Tevatron due to their large 
production cross sections at y/s — 1.96 TeV and the 
high integrated luminosity data sets available from the 
CDF and DO experiments during Run II. Detailed mea- 
surement of the W and Z properties at the Tevatron is 
not only important to further test GSW theory and the 
EWSB mechanism in the SM but also to search for new 
physics beyond the SM using the highest energy collisions 
currently available. We summarize the current Tevatron 
measurements of W and Z properties in Section IIII.AI 
through UlLD] 

The production of heavy vector boson pairs {WW, 
WZ, and ZZ) is far less common than inclusive W and 
Z production. While a W boson is produced in every 
3 million pp collisions and a Z boson in every 10 mil- 
lion, the production of a WW pair is a once in 6 billion 
event, WZ a once in 20 billion event, and ZZ a once 
in 60 billion event! Diboson production is sensitive to 
the triple gauge couplings (TGCs) between the bosons 
themselves via an intermediate virtual boson. The boson 
TGCs are an important consequence of the non-Abelian 
nature of the SM electroweak gauge symmetry group. At 
the highest accessible energies available at the Fermilab 
Tevatron, diboson production provides a sensitive probe 
of new physics, including anomalous trilinear gauge cou- 
plings, new resonances such as the Higgs boson, and large 
extra dimensions [12]. Recent results on diboson produc- 
tion from the Tevatron are discussed in Section IIII.EI 



A. Heavy Gauge Boson Production 

In high energy proton-antiproton collisions at the Fer- 
milab Tevatron, heavy vector bosons {V = W, Z) are 
produced at tree-level via quark-antiquark annihilation 
(qq — V) as shown in Fig. [TJ At high transverse momen- 




FIG. 1 Leading order (a) W and (b) Z boson production via 
quark-antiquark annihilation. The production is dominated 
hy q = u,d valence quarks. The 7* contribution interferes 
with the Z diagram but is not shown. 

tum^ , the leading-order QCD subprocesses are qq — >■ Vg 
and qg — >■ Vq. The production properties of heavy gauge 
bosons provide tests of perturbative QCD and, under cer- 
tain circumstances, information about quark and gluon 
momentum distributions within the proton and antipro- 
ton. 

The heavy vector boson production cross section in a 
pp collision is given by 

a{V) = j(jQyf{Xq)f{Xq)dXqdXq 

in which cro.i/ is the cross section for production of a vec- 
tor boson by a quark, antiquark pair with Feynman x 
values of Xq and x-q respectively, and f(xq) and f{x-q) 
are the parton distribution functions for the proton and 
antiproton. As defined here, contributions from pure 7* 
and Z/7* interference terms are not included but are ac- 
counted for in comparisons with theory in the measure- 
ments we describe. 

The remainder of this section includes a summary of 
current W and Z cross section measurements that the 
Tevatron and their comparisons with theory. As previ- 
ously mentioned, these measurements of W and Z boson 
production cross sections primarily a test QCD rather 
than electroweak theory. We include them here since de- 
tailed study of heavy gauge boson states are central to 
tests of electroweak theory and, therefore, it is impor- 
tant consider how well their production in pp collisions is 
understood. 

The cross section times branching fraction of W and Z 
bosons is measured in the fully-leptonic decay channels 
W ^ iv and Z — >• i£, where i = e, ^, r. While the final 
states involving r-leptons are important for many rea- 
sons, we restrict ourselves to the e and /i final states for 
the cross section discussion since these give the highest 



Throughout this paper, "transverse" is taken to be in a plane per- 
pendicular to the beam directions, and unless otherwise noted, 
quantities with a "T" subscript, e.g Et, are values projected onto 
this plane. Also, charge conjugation is assumed throughout. 



5 



precision measurements. The fully-leptonic decay chan- 
nels are chosen over hadronic channels for these mea- 
surements, since the latter suffer from large backgrounds 
to due to the hadronic decay of jets produced by QCD 
processes. 



1. W and Z Cross Sections 

Using a next-to-next to leading order (NNLO) predic- 
tion [III calculated with the MRST2004 NNLO parton 
distribution function [14] and the SM branching fractions 
for the W and Z bosons into fully-leptonic final states, 
the cross-section times branching fractions are calculated 
to be 

a{pp ~^W)x B{W tv) = 2687 ± 54 pb. 



and 



a{pp Z)x B{Z ££) = 251.9l^i°8 pb. 



The uncertainties are a combination of the MRST un- 
certainties and the difference between the central value 
above and that computed using the CTEQ6.1M parton 
distribution functions [15i |. 

The W and Z cross sections times branching fractions 
to fully-leptonic final states have been measured by CDF 
[l6l | using J C dt = 72 ph"^ . The results for electron and 
muon channels combined are 

a{pp ^W)x B{W -> Iv) = 

2749 ± 10(stat.) ± 53(syst.) ± 165(lumi.) pb 

and 

a{pp ^ Z)x B{Z ^ U) = 

254.9 ± 3.3(stat.) ± 4.6(syst.) ± 15.2(lumi.) pb. 

The measurements of the Z cross section include ad- 
ditional contributions from 7* and interference 
which give events that are experimentally indistinguish- 
able from the Z process. The size of these contri- 
butions depends on the mass range considered. 
For the range 60 GeY/c < M(Z/7*) < 130 GeV/c^ 
these contributions increase the cross section by a factor 
1.019 ±0.001 relative to the Z-only cross section, and for 
the mass range 66 GeV/c^ < M(Z/7*) < 116 GeV/c^, 
the cross section is increased by a factor 1.004 ± 0.001. 

A precision measurement of the ratio R of W and Z 
cross section times branching fraction given by 



R = 



a{pp -^W)x B{W if) 
aijip Z)x B{Z ei) 



can be used to test the SM. For example, new high mass 
resonances decaying to either W 01 Z bosons could lead 
to a deviation of the measured value of R from the SM 
expectation. 



Important systematic uncertainties such as the inte- 
grated luminosity uncertainty cancel in the measurement 
of R. The ratio R has been measured by CDF [1^ using 
J £ dt = 72 pb-i to be 

R = 10.84 ± 0.15(stat.) ± 0.14(syst.) 

This measurement has a precision of 1.9% and is con- 
sistent with SM expectation at NNLO of 10.69 ± 0.08 

El. 

A summary of the results are shown in Tab. HI The 
individual results are in good agreement with each other 
and with the prediction from theory. 

Even with the moderate sample sizes of these measure- 
ments, the W and Z cross section results are limited by 
the systematic uncertainty from the luminosity measure- 
ments. Because of this, further improvement in these 
measurements is not anticipated. CDF and DO Collab- 
orations continue to use Z-boson production to measure 
experimental efficiencies acceptances and cross-checks on 
temporal or instantaneous luminosity dependence of de- 
tector response. The updated cross section results could 
be reinterpreted as measurements of the integrated lumi- 
nosity. 



B. W Boson Mass 

At tree level, the W boson mass is fully-determined 
by the electromagnetic fine structure constant, the weak 
Fermi coupling and the cosine of the weak mixing angle. 
When higher-order EW corrections, like those shown in 
Fig. [21 are included, the expression is modified to [TtI 



naEM 



1 



V2Gf (1-M2,/M|)(l-Ar) 



(3) 



in which Mw (and later Tw) correspond to the parame- 
ters of a Breit-Wigner distribution with an s-dependent 
width. The term Ar includes the effects of radiative cor- 
rections and depends on and log Mh, where Mt and 
Mh are the top quark and Higgs boson masses, respec- 
tively. Measurement of the W boson mass can be used to 
constrain the allowed Higgs boson mass. The precisions 
of the W boson and top quark mass measurements are 
currently the limiting factors in the indirect constraint on 
the Higgs boson mass and are shown in Fig. [31 Earlier 
measurements of the W boson mass have been made by 
the LEP experiments and CDF ^ and DO ^ 

in Run I of the Tevatron. 

Signal-to-background and resolution considerations 
dictate use of the W ^ eu and W ^mv decays modes for 
the W boson mass measurement at the Tevatron. The 
momentum component of the neutrino along the beam 
direction cannot be inferred in pp collision events, so the 
W invariant mass cannot be reconstructed from its de- 
cay products and other variables are used to determine 
the mass. Three variables are used: (1) the lepton (e or 
/i) transverse momentum p^, (2) the (inferred) neutrino 
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Integrated 
Luminosity 


Data 
Yield 


Predicted 
Background 


Axe 


Measured (j x Br (pb) 


W -^eiy 
W ^ fj-v 
and fiiy 


72 pb-i 
72 pb^^ 


37584 
31722 


1762 ± 300 
3469 ± 151 


U.± ( »^_o.0038 
U.i^'iZ.Q 0034 


2.771 ± O.OUto^otl ± 0.166 nb 
2.722 ± 0.015t;];E;^? ± 0.163 nb 

7/1 Q -t-A nin -1-0 nf^Q -i-n nK 
z. ( 4y it u.uiu zt u.uoo it u. loo iid 


Z ^ ee 
Z ^ Hfi 
ee and ^fj, 


72 pb-^ 
72 pb"^ 


4242 
1785 


62 ± 18 
13± 13 


U.ZZDy_o.o048 

n992+° °''28 
u.uyyz_o 0031 


255.8 ±3.9t5:4± 15.3 pb 
248.0 ± 5.9l?-^ ± 14.8 pb 
254.9 ± 3.3 ± 4.6 ± 15.2 pb 



TABLE I The measured integrated luminosity, yield, background, signal acceptance times efficiency and resulting cr x Br for 
each of the four Z cross section times branching ratio analyses by the CDF Collaboration in 46]. The combination of the 
dielectron and dimuon channels using the BLUE method is also shown. The Z/"/* results have been corrected to Z only as 
described in the text. 





W / W 

VWWWNTWVAAAAAAA/V/V/ 

FIG. 2 Feynman diagrams showing corrections to the W bo- 
son mass from (a) the top quark and (b) the SM Higgs boson. 



transverse momentum and (3) the transverse mass, 
rriT = \J p^j.p^[l — cos((5(^)] in which 5(j) is the lepton- 
neutrino opening angle in the plane transverse to the 
beam. Although these variables are highly correlated, 
their systematic uncertainties are dominated by different 
sources, so the results are combined taking their statisti- 
cal and systematic correlations into account. 

The event decay kinematics in the transverse plane 
are fully characterized by two quantities: (1) the lepton 
transverse momentum and (2) the transverse momen- 
tum of the hadronic recoil ut required to balance the 
transverse momentum of the the W . The hadronic recoil 
is defined as the vector sum of all energy measured in 
the calorimeter excluding that deposited by the lepton. 
From these two measurements, the neutrino transverse 
momentum is inferred: —p^ = p^ + iir. In practice, 
the recoil is computed using the event transverse missing 
energy (.^y) measured in the calorimeter after removing 
the contribution to calorimeter energy associated with 
the lepton. The W boson mass is determined by gen- 
erating predicted distributions (templates) of the three 
measurement variables for a range of input W boson 
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m^ [GeV/c2] 



FIG. 3 The measured top quark and W boson masses and 
a band of allowed Higgs boson masses. [2^ This includes the 
recent results summarized in the text. 



mass hypotheses. These are generated using dedicated 
fast Monte Carlo simulation programs, and the mass is 
determined by performing a binned maximum likelihood 
fit of these templates to the distributions observed in 
data. The Run II measurements use a blinding tech- 
nique in which an unknown offset is added to the value 
returned from the fits until the analyses are finished. At 
that point, the offset is removed to reveal the true value. 
The results from the different fit variables {p^, p^, m-r) 
and final states (e or p) are combined using the BLUE 
algorithm [Ullil. 

Both DO and CDF have reported mass measurements 
using Run II data. The CDF result fl!] uses J £ dt = 
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CDF 



W 



fC dt 
E^, cal 
Pt, trk 

M 

Ut 

Yield 

Background 



0.2 fb"^ 
> 30 GeV 

> 18 GeV/c 

> 30 GeV/c 
< 15 GeV/c 

63964 
7.5% 



W 
0.2 fb^ 



> 30 GeV/c 

> 30 GeV/c 
< 15 GeV/c 

51128 
1.1% 



DO 
W ^ eiy 
1.0 fb"^ 
> 25 GeV 

> 10 GeV/c 
< 1.05 

> 25 GeV/c 
< 15 GeV/c 

499830 
4.0% 



TABLE II The integrated luminosity, kinematic selection, 
event yield and backgrounds for the CDF and DO W boson 
mass analyses. The background is given as a percentage of 
the total yield. 



0.2 fb ^ and results are reported for W eiy and W ^■ 
fiiy decays. The DO result [13 uses J C dt ^ 1.0 fb^i 
and results are reported for the W ^ ev decay mode 
only. Candidate events are selected by requiring a single 
high-pT isolated charged lepton and large Tab. 
shows the kinematic selection requirements, event yields 
and background fraction for the W event selections. 

The backgrounds include Z ^ il events in which one 
lepton escapes identification, WW diboson events, and 
misidentification backgrounds in which the lepton is ei- 
ther a jet misidentified as an electron or a muon from 
semileptonic decay of hadrons in which the rest of the 
associated hadronic jet is not reconstructed. An addi- 
tional source of events are the sequential decays W — ?> 
TV — >■ evvv and W ^ tv ^ [ivvv. The CDF analysis 
treats these as signal while the DO analysis incorporates 
these into the background template distributions. 

The in situ calibration of charged particle momenta 
(CDF) and calorimetric measurement of electron energy 
(CDF, DO) is of crucial importance to this result. The 
CDF analysis uses a calibration of the tracker momen- 
tum scale {p) determined from dimuon and dielectron 
decays of J/ ip, T and Z particles. This calibration is then 
transferred to the calorimeter energy measurement {E) 
using the E/p ratio. A final improvement is made for the 
W ^ eu mode by incorporating an additional calorime- 
ter calibration based on Z — )■ ee decays. The DO analysis 
uses calorimeter energy measurements, and the calibra- 
tion is based on the mass reconstructed in Z — >■ ee events 
and a detailed simulation of the calorimeter response. 
For both experiments, this calibration is the dominant 
source of systematic uncertainty. Other sources of sys- 
tematic uncertainty arise from trigger efficiency, lepton 
identification efficiency, correlation (in)efficiency such as 
occurs when the hadronic recoil is near the charged lep- 
ton, backgrounds, electroweak and strong contributions 
to the production and decay model and the parton dis- 
tribution functions. 

The rtiT distributions for each channel are shown in 
Fig. UJ and the results from each channel and the combi- 
nations of the channels for each experiment are shown in 
Tab. mil The systematic uncertainties are dominated by 
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FIG. 4 The mr distributions for (a) the CDF W ^ chan- 
nel, (b) the CDF W ^ ev channel and (c) the DO W eu 
channel. 



the lepton energy calibration. This contributes 17 MeV 
uncertainty to the CDF and p^ channels, 30 MeV to 
the CDF rriT channel and 34 MeV to each DO channel. 

These results are shown along with previous measure- 
ments in Fig. [SJ The world average combination [24; 2^ 
has been updated with these measurements using the 
BLUE [2^ [2^ algorithm including correlations. The re- 
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Mw (MeV) 


xVdof 


CDF W -¥eu 


rriT 


80493 ± 48 ± 39 


86/48 




Pt 


80451 ± 58 ± 45 


63/62 




$T 


80473 ± 57 ± 54 


63/62 


Combined 




80477 ± 62 




CDF W ^ iM/ 


rriT 


80349 ± 54 ± 27 


59/48 




M 

Pt 


80321 ± DO ± 4U 


72/62 




jp 

ipT 


8U39D ± DO ± 4d 


/I A /GO 

44/03 


Combined 




80352 ± DO 




CDF W ^eu + fiu 




80413 ± 34 ± 34 




mW ^ev 


rriT 


80401 ± 23 ± 37 


48/49 




Pt 


80400 ± 27 ± 40 


39/31 




$T 


80402 ± 23 ± 43 


32/31 


DO Combined 




80401 ± 21 ± 38 





TABLE III The individual CDF and DO W boson mass re- 
sults and their combinations. When two uncertainties are 
given, the first is the statistical uncertainty and the second is 
the systematic uncertainty. 



CDF Run 0/1 
DO Run I 
CDF Run II 
DO Run II 

Tevatron 2009 

LEP2 average 
World average 



80.436 ± 0.081 [24] 
80.478 ± 0.083 [25] 
80.413 ±0.048 [28] 
80.402 ± 0.043 [29] 

80.420 ± 0.031 [30] 

80.376 ± 0.033 [20-23] 
80.399 ± 0.023 [30] 



L 



80 80.2 80.4 

m^^ (Gey Id) 



80.6 



FIG. 5 Ideogram of previous measurements, the Run II CDF 
and DO measurements and the new world average. The Ref. 
is shown to the right of each measurement as [XX] . 



C. W Width 

Like the W boson mass, the width Tw is also predicted 
by the SM. It is given by 



Tw = (3 + ifqcn) 



6y27r 



(l + <5) 



(4) 



in which focD ~ 3.1 is a QCD correction factor, and 
6 — 2.1%^] is an EW correction factor. Direct measure- 
ments of Tw were made by CDF 3C | and DO [3l| using 
Tevatron Run I data and combined [31|. Measurements 
were also made by the experiments at LEP [l8- 21]. 

Because the W boson mass is distributed according to 
a Breit-Wigner, there is a tail at large mass values. The 
W boson width result is obtained using the rriT distri- 
bution in a region where the shape and event yield are 
dominated by events from the high mass region of the 
Breit-Wigner with limited impact from detector resolu- 
tion effects. As for the W boson mass measurement, a 
binned likelihood comparison of the observed spectrum 
to templates generated for different W boson widths is 
used to extract the numerical value. 

The W width was measured by CDF [H and DO [H 
using Run II data. The CDF result uses a data set of 
J £ dt — 0.35 fb~^ and the electron and muon final 
states. The DO result uses J £ dt — 1 fb~^ and the 
electron final state. The DO data set is the same one 
used for the W boson mass measurement. Similar ded- 
icated simulations and processing were used for the W 
width measurement as were used for the W boson mass 
measurement. DO used a different hadronic recoil pro- 
cedure [33 |. The new procedure gives a W boson mass 
result consistent with the standard method. Fig. [5] shows 
the high mass region of the rriT distribution for the three 
channels analyzed: the CDF W ^ fjiv channel, the CDF 
W ^ ev channel and the DO W ^ ev channel. The Tw 
results are given in Tab. IIVI and shown in Fig. [71 The 
systematic uncertainties are dominated by hadronic re- 
coil scale and resolution uncertainties. These contribute 
54(49) MeV for the CDF W ^ ev {W ^ ^lu) channel 
and 41 MeV to the DO result. Other important sources 
include the lepton scale uncertainty and background un- 
certainty which are one half to two thirds the size of the 
hadronic recoil uncertainty. 



suit is 

Mw = 80.420 ± 0.031 GeY/c. 

Because the systematic uncertainties are dominated by 
the statistical precision of the calibrations determined 
from control data samples, the systematic uncertainty in 
future measurements is expected to improve as the inte- 
grated luminosity increases. The ultimate limiting sys- 
tematic is expected to be that introduced by the parton 
distribution functions. In the current results, this con- 
tributes an uncertainty of 11 MeV to all channels with a 
100% correlation among the channels. 



D. Forward-backward Asymmetry, Afb 

Production of Z bosons at the Tevatron is dominated 
by the process qq — > Z/^* in which q — u,d are pro- 
ton valence quarks. The SM couplings of the Z and 7 
to quarks depend on the quark charge Q, the isospin /a, 
and the sine of the weak mixing angle sin^ 6*14^- The dif- 
ferential cross section as a function of the direction of the 
fermion resulting from the Z/^* decay is given by 

, ^ = a(l + cos^ 6) + b cos 9 (5) 
a cosy 
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Channel 


Yield 


Fit Range(GeV/c^) 


Tw (MeV) 


xVndof 


CDF W ^ ^if 


2619 


90 < -m-T < 200 


1948 ± 67 ± 71 


17/21 


CDF W -^eiy 


3436 


90 < mr < 200 


2118 ± 60 ± 79 


19/21 


CDF Combined 






2032 ± 45 ± 57 






5272 


100 < mr < 200 


2028 ± 39 ± 61 


75.2/75 



TABLE IV The Tw measurements. For the result, the first uncertainty is the statistical uncertainty and the second is the 
systematic uncertainty. 



in which 9 is the angle of the fermion from the Z/j* decay 
measured relative to the incoming quark direction in the 
Z/"f* rest frame. The relative Z and 7* contributions to 
the cross section vary as a function of the Z/7* mass, 
and differences in the Z and 7 couplings to quarks result 
in different angular distributions for the decay products 
for up-type {I3 = +1/2) and down-type {I3 — ^1/2) 
quarks. Together, these two effects produce mass and 
flavor dependence in the coefficients a and b which can 
be calculated assuming the SM. 

The forward-backward Z/j* production charge asym- 
metry is defined as 



Dielectron Mass 



A = 



(j^ — (J- 



in which (t+ and (7_ are the integrated cross sections 
for the cases cos 61 > and cos6' < respectively. The 
asymmetry extracted experimentally is given by 



A 



FB 



N+ - N. 
N+ + N. 



(6) 



in which Nj^ and N- are the acceptance, efficiency 
and background corrected fermion yields in the forward 
(cos0 > 0) and backward {cosO < 0) directions respec- 
tively. Measuring the asymmetry rather than differen- 
tial cross sections allows cancellation of many systematic 
uncertainties, particularly those affecting the overall nor- 
malization. 

Measurements of the asymmetry as a function of di- 
electron mass have been made by both CDF [s^ and 
DO [3^ using the dielectron final state. The CDF result 
uses a sample with J C dt = 72 pb~^, and the DO mea- 
surement uses J C dt = lA fb^^. The selection criteria 
are similar to those for the Z cross section measurements 
although a larger dielectron mass range was selected for 
the asymmetry measurements. Two experimental issues 
of particular importance to this measurement are (1) con- 
trolling asymmetries in either detector acceptance or se- 
lection efficiency as a function of dielectron mass and (2) 
limiting the impact of electron charge misidentification. 

Tab.fVland Tab. I VII give the dielectron mass range and 
the predicted and measured values of ApB for each mass 
bin, and Figs. [8] and [9] show the measured asymmetries 
and the SM predictions as a function of mass for the CDF 
and DO results respectively. 

Using these measurements and the SM prediction for 
the coefficients a and 6, sin^ ^^'^ can be determined. 



ige 


(GeV/c2) 


Pythia 


Measured 


40 


-60 


-0.214 ±0.003 


0.11 ±0.13 ±0.05 


60 


- 70 


-0.420 ± 0.005 


0.5ltoi? ±0.07 


70 


- 78 


-0.410 ± 0.005 


0.45 ±0.19 ±0.05 


78 


-86 


-0.214 ±0.003 


0.11 ±0.17 ±0.09 


86 


-88 


-0.079 ± 0.004 


0.07 ±0.23 ±0.19 


88 


-90 


-0.001 ±0.002 


0.03 ±0.13 ±0.09 


90 


-92 


0.054 ±0.001 


0.047 ± 0.077 ± 0.076 


92 


-94 


0.112 ±0.002 


0.15 ±0.19 ±0.16 


94 


- 100 


0.198 ±0.003 


0.35 ±0.20 ±0.08 


100 


- 105 


0.338 ± 0.006 


0.02 ±0.30 ±0.19 


105 


- 120 


0.454 ± 0.006 


0.67 ±0.15 ±0.05 


120 


- 140 


0.554 ± 0.002 


0-S2to ll ±0.01 


140 


-200 


0.598 ± 0.002 


0.29tg:i? ±0.01 


200 


-300 


0.609 ± 0.004 


0.65+°oi° ± 0.03 


300 


-600 


0.616 ±0.007 


i.uu_Q 5i_Q 02 



TABLE V The expected and measured asymmetries as a 
function of dielectron mass (CDF). 



Dielectron Mass 






Afb 


Range 


(GeV/c)2 


Pythia 


ZGrad 


Measured 


50 


-60 


0.293 


0.307 


0.262 ± 0.066 ± 0.072 


60 


- 70 


0.426 


0.431 


0.434 ± 0.039 ± 0.040 


70 


- 75 


0.449 


0.452 


0.386 ±0.032 ±0.031 


75 


-81 


0.354 


0.354 


0.342 ± 0.022 ± 0.022 


81 - 


-86.5 


0.174 


0.166 


0.176 ±0.012 ±0.014 


86.5 


- 89.5 


0.033 


0.031 


0.034 ± 0.007 ± 0.008 


89.5 - 92 


0.051 


0.052 


0.048 ± 0.006 ± 0.005 


92 


-97 


0.127 


0.129 


0.122 ±0.006 ±0.007 


97 


- 105 


0.289 


0.296 


0.301 ±0.013 ±0.015 


105 


- 115 


0.427 


0.429 


0.416 ± 0.030 ± 0.022 


115 


- 130 


0.526 


0.530 


0.543 ± 0.039 ± 0.028 


130 


- 180 


0.593 


0.603 


0.617 ±0.046 ±0.013 


180 


- 250 


0.613 


0.600 


0.594 ±0.085 ±0.016 


250 


- 500 


0.616 


0.615 


0.320 ±0.150 ±0.018 



TABLE VI The expected and measured asymmetries as a 
function of dielectron mass (DO). For the measured values, 
the first uncertainty is statistical and the second is systematic. 



Here 0^'^ is the weak mixing angle including higher order 
corrections. The current world average is 



qeff 



0.23149 ±0.00013 



using the AIS scheme [37|. Among the measurements 
used for the world average are two, the charge asymme- 
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r„ = 1948±67MeV 
X^/dof [fit range] =17/21 
X'/dof [full range] = 21/29 




Data 

W ^ [J.V MC + Background 

|»j Background: 5.7%z^^^x 

' 2.0% W -> TV 

0.3% Multijet 
0.2% DIF 



150 200 
m-r ()xv)(GeV/c') 



Standard Model 



DO Runll (e) 



2.028 ± 0.072 [34] 



CDF Runll (e+n) 2.032 ± 0.073 [33] 

DO RunI (e) 2.231 ± 0.173 [32] 

CDFRunl(e+n) 2.05010.128 [31] 

Had. Col. Avg (2008) 2.056 ± 0.062 [32] 



Prel. LEP Avg (2006) 2.1 96 ± 0.084 [20-23] 
World Avg (2008) 2.1 06 ± 0.050 [38] 



1.8 2 2.2 2.4 2.6 2.8 3 
Direct Measurements (GeV) 



3.2 




m-r (ev)(GeV/c') 



FIG. 7 Summary of the Tw measurements. The Hadron Col- 
hder Average and World Average results do not include the 
DO [H result. 
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150 200 
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FIG. 6 The rriT distributions used for the Fvi^ measurement, 
(a) The CDF W ^ fiu channel, (b) The CDF W ^ eiy 
channel, (c) The DO W ev channel. 



FIG. 8 Dielectron forward-backward asymmetry as a function 
of dielectron mass(CDF). 



try for 5-quarlc production |Tl| from LEP and SLD and 
the measurement from NuTeV [s^ , which differ from the 
world average by more than two standard deviations. 

The values for sin^ 6*^-^ extracted using fits to the CDF 
and DO ApB distributions are 

sin^ 0^^ = 0.2238 ± 0.0040 (stat) ± 0.0030 (syst) 

for CDF and 

sm^0y(r^ = 0.2327 ± 0.0018 (stat) ± 0.0006 (syst) 

for DO. Fig. [To] shows these results compared to other 
measurements. The results from DO are comparable in 
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FIG. 9 Dielectron forward-backward asymmetry as a function 
of dielectron mass (DO). 
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FIG. 10 Comparison of the Tevatron asymmetry results with 
those from other experiments. 
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FIG. 11 The Zuu and Zdd coupling constants. [35 



precision to other measurements for light quarks. The 
current Tevatron results are limited by sample statistics, 
but by the end of the Tevatron running, CDF and DO 
are expected to have the most precise measurements of 
sin^ 9ll^ for light quarks. 

CDF also removed the assumption of SM quark cou- 
plings and determined values from a four parameter fit of 
ApB measurements to the SM prediction as function of 
the vector and axial vector couplings for u and d quarks. 
The fit has a x^/dof = 10.4/11, and the resulting cou- 
pling values are shown (with the SM values) in Fig. [11] 
No evidence of deviation from the SM is observed. 



E. Dibosons 

1. Trilinear Gauge Couplings (TGCs) 

The non-Abelian nature of the gauge theory describing 
the electroweak interactions leads to a striking feature of 
the theory. In quantum electrodynamics, the photons 
carry no electric charge and thus lack photon-to-photon 
couplings and do not self-interact. In contrast, the weak 
vector bosons carry weak charge and do interact amongst 
themselves through trilinear and quartic gauge boson ver- 
tices. Fig. [T2| shows the tree-level diagram for diboson 
production involving trilinear gauge couplings. 
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7, 2", 11 



FIG. 12 Leading order diagram for diboson production via 
quark-antiquark annihilation involving the trilinear gauge 
coupling. 

The SM Lagrangian that describes the WWV{V = 
Z, 7) interaction is given by 

where W^' denotes the W field, W^^^ = d^^W,, - d^W^, 
V^u = dyY^ — d^V^, the overall couplings are gwwt = — e 
and gwwz = — ecot^vi'j and is the weak mixing 
angle [3^. At tree level in the SM, the trilinear boson 
couplings involving only neutral gauge bosons (7 and Z) 
vanish because neither the photon nor the Z boson carry 
electric charge or weak hypercharge. 

A common approach used to parameterize the low en- 
ergy effects from high-scale new physics is the effective 
Lagrangian approach that involves additional terms not 
present in the SM Lagrangian ^39] . This approach is con- 
venient because it allows for diboson production proper- 
ties measured in experiments to be interpreted as model- 
independent constraints on anomalous coupling parame- 
ters which can be compared with the predictions of new 
physics models. 

A general form for the WWV Lorentz-invariant inter- 
action Lagrangian with anomalous coupling parameters 
gY , KVi and Ay is given by [sol: |40| 

C'^wv = ^gwwv[gX{w+.w-^ -W+^^W;,)V'' 

Note that reduces to C^^ for the values A..^ = Az = 
and gj — gi — — Kz = 1. Deviations from the 
SM values of the coupling parameters are denoted by 
AgY , AKy, and AAy. We have assumed that C and P 
are conserved in the interaction Lagrangian. There is no 
reason to believe that this assumption is valid unless the 
physics that leads to anomalous couplings respects these 
symmetries. It is straightforward to include additional 
terms that violate C and P, but we refrain from doing 
so in order to keep the discussion simple. 

Electromagnetic gauge invariance requires AgJ = 0. 
The W boson magnetic moment fiw and the electric 



quadrupole moment Qw are related to the coupling pa- 
rameters by 

fiw = 7i (1 + K7 + A.^) 

and 

Qw — 9~('*7 ~ -^7) 

The anomalous couplings (aside from gj) are usually 
assumed to have some dependence on an energy scale 
(form factors) which suppresses them at large scales to 
avoid violation of tree-level unitarity in the diboson pro- 
duction amplitude (4TI : . The parameterization gener- 
ally used for the energy dependence of a given coupling 
parameter a is 

^ ' (l + .s/A2)2 

where Vs is the partonic center-of-mass collision energy, 
ao is the value of the coupling parameter in the limit 
s — > 0, and A is the cutoff scale. When deriving cou- 
pling limits from hadron collider data, the value of A is 
taken to be close to the hadron collision energy; even 
large variations (e.g. 50%) of A around this scale have 
minimal impact on the results. Physically, the scale A 
can be considered the scale at which the new physics re- 
sponsible for the anomalous coupling is directly accessible 
(e.g. through pair production of new particles) . This ap- 
proach is different from the effective field theory approach 
discussed in Section |lll where the coefficients of higher- 
dimension operators are constants. While in the same 
spirit as effective field theory, the effective Lagrangian 
approach to anomalous couplings is different in practice. 
In particular, the form factors invoked in the effective 
Lagrangian approach are unnecessary in an effective field 
theory approach. 

In the presence of new physics, neutral TGCs (those 
involving only 7 and Z bosons) can contribute to Zj and 
ZZ production. As previously described, neutral TGCs 
are anomalous by their very nature since these couplings 
are absent in the SM. For each of the diboson final states 
Z7 and ZZ, one can follow an analogous procedure to the 
anomalous charged TGCs (those involing a W boson) by 
writing down the most general effective Lagrangian that 
respects Lorentz invariance and electromagnetic gauge 
invariance [H; H^]. Using prescriptions detailed in [4^ 
and d^l, the effective Lagrangians introduce anomalous 
coupling parameters hfoiV = ^,Z and i = 3,4) and 
fj^{V = j,Z and j = 4,5) respectively, which can be 
constrained through an analysis of Zj and ZZ produc- 
tion in high-energy collider data. It is important to note 
that, under the assumption of on-shell Z bosons, the Z"fZ 
couplings contributing to Zj production and ZZ7 cou- 
plings contributing to ZZ production are completely in- 
dependent |44]. 

In general, the effect on observables from turning on 
anomalous TGCs are correlated. When we refer to "ID 
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limits," we refer to the limits derived on one parameter 
when the others are set to their SM values. 

There are a few important differences regarding the 
study of diboson physics in particle collisions at LEP, 
Tevatron, and the LHC that are worth pointing out at 
this stage: 

• At LEP, e+e~ collisions occur at a well-defined en- 
ergy that is set by the accelerator. Therefore, the 
center-of-mass energy is known with good precision 
and there are no form factors in anomalous coupling 
analyses. 

• In e+e^ collisions, the initial state has zero elec- 
tric charge. Therefore, exclusive states with net 
charge, such as WZ and W^, cannot be produced 
at LEP. The WW and ZZ states can and have been 
produced and studied at LEP. A measurement of 
the WW cross section over a scan in beam energy 
dramatically illustrates the existence of the WW Z 
coupling in electroweak theory [45j . 

• At hadron colliders, ^/s is fixed for long periods of 
time (defining different periods of the accelerator 
operation that change very infrequently) but \/l 
varies collision by collision. Any anomalous cou- 
plings are likely to be a/I dependent. The form 
factor ansatz used to cut off the anomalous cou- 
pling parameters at large a/I to preserve S-matrix 
unitary of the amplitude reflects this dependence. 
For this reason, it is reasonable to expect that, 
once a sufficient amount of integrated luminosity 
has been acquired, the higher-energy reach afforded 
by high-energy hadron collisions will lead to better 
sensitivity to anomalous couplings as compared to 
the limits from LEP, despite larger backgrounds in 
a typical hadron collision event. In other words, 
hadron collisions at the Tevatron and LHC sam- 
ple events with a larger average a/I as compared 
to LEP collision energies and it is exactly those 
high \f% events that are most sensitive to effects 
of anomalous couplings from new physics at higher 
energy scale. 

• Because the Tevatron is a collider, the pro- 
duction cross sections for W~^Z and W^Z are 
equal. The same is true for and W~^ pro- 
duction. When produced in a pp collider such as 
the LHC, positive and negative net charge diboson 
states have different production cross sections (e.g. 
a{pp W+Z) > a{pp W'Z)). 

2. W7 

The W^ final state observed at hadron colliders pro- 
vides a direct test of the WW"/ TGC. Anomalous WWj 
couplings lead to an enhancement in the production 
cross section and an excess of large Et photons. Both 





DO Analysis 
ez^7 ^wy 


CDF Analysis 


J£dt (fb-^) 


0.16 


0.13 


20 (0 17) 


19 fO 181 


W+ jet(s) 


59 ±5 


62 ±5 


60 ±18 


28 ±8 


ieX 


1.7±0.5 


0.7±0.2 








0.42 ±0.02 


1.9±0.2 


1.5 ±0.2 


2.3 ±0.2 






6.9 ±0.7 


6.3 ±0.3 


17.4 ±1.0 


Total Bkg. 


61 ±5 


71 ±5 


67 ±18 


47 ±8 


^ observed 


112 


161 


195 


128 


o-x BR (pb) 


13.9 ±3.4 


15.2±2.5 


19.4 ±3.6 


16.3 ±2.9 


Theory ctx BR 


16.0 ± 0.4 


19.3 


± 1.4 



TABLE VII Summary of the W7 ^ tv-i CDF [H and DO 
cross section analyses compared with the theory [43| ) . For 
the CDF analysis, the luminosities with and without paren- 
theses correspond to central and forward leptons, respectively. 
The theory cross section is the NLO calculation (in pb) from 

CDF and DO have published measurements of the W^ 
cross section using leptonic decays of the W bosons and 
J Cdt = 0.2 fb"^ [Slii- The signature of the W-f signal 
is an isolated high Et lepton, an isolated high Et photon, 
and large $!t "iue to the neutrino from the W decay. The 
dominant background is from M^+jets where a jet mimics 
an isolated photon. A lepton-photon separation require- 
ment in 77 — space of AR = ^ (ATy)^ + (A^)^ > 0.7 
is made by both CDF and DO to suppress events with 
final-state radiation of the photon from the outgoing lep- 
ton and to avoid coUinear singularities in theoretical cal- 
culations. A kinematic requirement on photon Et of 
Et > 7(8) GeV is made by CDF (DO ) in the analysis. 
CDF measures 

a{pp -^Wj + X)x BR{W Iv) 

^ 18.1 ± 1.6(stat.) ± 2.4(syst.) ± 1.2(lum.) pb 

[4^ in agreement with the next-to- lea ding order (NLO) 
theoretical expectation [Et > 7 GeV) [isf of 19.3 ± 1.4 
pb. DO measures 

a{pp ^W"/ + X)x BR{W Iv) 

= 14.8 ± 1.6(stat.) ± 1.0(syst.) ± 1.0(lum.) pb 

[47I I also in agreement with the NLO expectation {Et > 8 
GeV) (43 of 16.0 ± 0.4 pb. Tab. EII] summarizes the 
Wj cross section measurement results. 

Both of the CDF and DO Wj cross section measure- 
ment are consistent with the SM expectations at NLO. In 
a more recent analysis .48] , DO uses four times more inte- 
grated luminosity as compared to [l^l and adds photons 
reconstructed in their endcap calorimeters (1.5 < \r]det\ < 
2.5) to search for anomalous WWj couplings based on 
the observed photon Et spectrum [i^l for photons with 
Et > 9 GeV. Additionally, the three-body transverse 
mass of the photon, lepton and J^t must exceed 120 
GcV/c (110 GeV/c) for the electron (muon) channel 
in order to suppress final state radiation. The photon 
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FIG. 13 The (a) photon Et distribution for the DO Wj anal- 
ysis The photon Et distribution of events are used to 
constrain anomalous couphng parameters shown in (b). 



Et spectrum and anomalous TGC limits are shown in 
Fig. [131 A LO simulation i49|] of the signal is used 
with NLO corrections [13] applied to the photon Et spec- 
trum. The one-dimensional limits at 95% confidence level 



(CL) are -0.51 < Ak^ < 0.51 and 



-0.12 < < 0.13 



for A = 2.0 TeV. 

The SM W^7 production involves interference between 
the amplitudes for a photon radiated off of an incoming 
quark (QED initial state radiation) and the photon pro- 
duced from the WWj vertex. This interference leads to 
a zero amplitude for the SM in the photon angular dis- 
tribution [5ll - l54l |. In W-f production, the radiation am- 
plitude zero (RAZ) manifests itself as a dip at ~ —1/3 in 
the charged-signed rapidity difference Qe x Ay between 
the observed photon and the charged lepton from decay 
of the W boson [s^ . Experimentally, the pseudorapidity 
difference Arj is used in place of the rapidity difference 
Ay, since it involves only the production angle with 
respect to the beam line (77 — — ln(tan(6'/2))) and is a 
very good approximation to Ay in the limit of massless 
particles. Using the same data they used to limits on 
WW J anomalous TGCs, DO made a first detailed study 
of the X Ay to search for the RAZ effect . 

Fig. [TiT a) shows the x Ay distribution of data com- 



pared with the SM expectation, which has a 15% prob- 
ability for compatibility between the data and the SM 
expectation, demonstrating a reasonable level of agree- 
ment. To specifically investigate the dip region around 
Qi X Ay = —1/3 indicative of the RAZ effect, a sim- 
ple test statistic R is constructed which is the ratio of 
the number events observed in a bin including the dip 
region to the number of events observed in a bin with 
more negative Qi x Ay where a maximum is expected, 
based on Monte Carlo, from the SM including the total 
DO acceptance in this analysis. A value of i? = 0.64 is 
observed in the data and it is determined that 28% of SM 
pseudo-experiments give a value as large or larger than 
0.64. 

A particular anomalous coupling hypothesis corre- 



sponding to 



0, 



-1 is chosen to test against 



as a null hypothesis since it is a specific model which 
leads to a Qe x Ay distribution without a dip as shown 
in the inset of Fig. [T4l b). Only the shape is tested such 
that the null "no-dip" distribution is normalized to the 
yield expected from the SM. The R distribution for the 
SM and no-dip hypotheses are shown in Fig. [T4l^b) along 
with the value (0.64) observed in the data. From this 
analysis, the no-dip hypothesis is excluded at 2.6cr Gaus- 
sian significance. 



3. Z-y 



As mentioned in Section IIII.E.ll photons do not di- 
rectly couple to Z bosons at tree-level in the SM. There- 
fore, observation of such a coupling would constitute evi- 
dence for new phenomena. The Z'y final state at hadron 
colliders involves a combination of ZZj and .^77 cou- 
plings. Both CDF and DO have made measurements of 
the Zj cross section in leptonic decay channels of the Z 
boson. The signature of the Zj signal is two isolated 
high Et charged leptons having the same flavor and op- 
posite charge with invariant mass consistent with decay 
of a Z boson, and an isolated high Et photon. The 
dominant background is from Z-|-jets where a jet mimics 
an isolated photon. As in the 11^7 analyses described in 
Section IIII.E.21 a lepton-photon separation requirement 
is imposed. A requirement of AR > 0.7 is made by both 
CDF and DO. Kinematic requirements on the photon Et 
of Et > 7(8)GeV/c^ and the dilepton invariant mass of 
Mu > 40(30) GeV/c^ are made by CDF (DO ) in the 
analysis. 

Using fCdt = 0.2 fb"\ CDF measures 

a{pp Z'j + X)x BR{Z II) 

= 4.6 ± 0.5(stat.) ± 0.2(syst.) ± 0.3(lum.) pb 

[4^ in agreement with the NLO theoretical expectation 
using Ref. [13] of 4.5 ± 0.3 pb using the CDF acceptance. 
Using / C dt — l.l fb^^, DO measures 

(j{pp Zj + X)x BR{Z -^■ II) 

= 5.0 ± 0.3(stat. -f syst.) ± 0.3(lum.) pb 
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FIG. 14 (a) The background-subtracteci charge-signed rapid- 
ity diflerence of events is shown in along with the SM ex- 
pectation from [isl ]. Consideration of the full covariance ma- 
trix gives a 15% probability for compatibility between the 
data and the SM expectation, (b) Distributions of the R-test 
statistic for the SM (solid line) and "no-dip" (dashed line) 
hypothesis pseudo-experiments from [i^. The data result is 
0.64 shown by the vertical line where only 45 out of 10* no- 
dip pseudo-experiment had a value smaller as small or smaller 
than 0.64. The R-test statistic is defined as the ratio of events 
in the "Dip Bin" to the "Small Max Bin" , as shown in the in- 
set where the solid line is the SM expectation and the dashed 
line is for the no-dip hypothesis. 



[56| also in agreement with the NLO expectation using 
the generator described in [57] of 4.7 ± 0.2 pb and the 
DO acceptance. Tab. I Villi summarizes the cross sec- 
tion measurement results. Fig. [15] shows the three-body 
mass Mu^ versus the dilepton mass Mu for candidate 
events in the DO analysis [56| . 

Using J C dt — 3.6 fb^^, the DO collaboration has ob- 
served (5.1cr) the Zj — > viy^ process for the first time 
at a hadron collider [5g| and used these events to search 
for Zjj and ZZ'y couplings that are absent at tree-level 
in the SM. The experimental signature for Zj — >■ lyi'j 
is a high energy photon and large $!rp. In the analysis, 
events are required to have a single photon candidate 
with Et > 90 GeV and > 70 GeV/c. The primary 



DO Analysis 



CDF Analysis 



J£dt (fb-^) 

Z+ jets background 

A^obscrvcd 

crx BR (pb) 
Theory crx BR 



1.1 1.0 

55 ± 8 61 ± 9 
453 515 
5.0 ± 0.4 
4.7 ± 0.2 



0.20 (0.17) 0.19 (0.18) 
2.8 ±0.9 2.1 ±0.6 

36 35 
4.8 ±0.9 4. 4 ±0.8 
4.5 ± 0.3 



TABLE VIII Summary of the Z-y CDF l46|Und DO ^ 

cross section analyses compared with the theory [53 • For the 
CDF analysis, the luminosities with and without parentheses 
correspond to central and forward leptons, respectively. The 
theory cross section is the NLO calculation (in pb) from [stI ]. 
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FIG. 15 The three-body mass Mu-y versus the dilepton mass 
Mii for Z'y candidate events in the DO analysis [5a | 



backgrounds are from W ^ ev and events unrelated to 
the collision in which muons from beam halo or cosmic 
rays produce energetic photons through bremsstrahlung. 
The W ^ eu is suppressed by removing events with 
high-pT tracks and the non-collision events are removing 
using available zq production information from the EM 
calorimeter and pre-shower detectors. A summary of the 
background estimates and observed events is shown in 
Tab. IIXI The measured cross section is 

a{pp Zj + X)x BR{Z w) 

= 32 ± 9(stat. ± syst.) ± 2(lum.) fb 

for the photon Et > 90 GeV, consistent with the NLO 
cross section of 39 ± 4 fb \5l\ . 

Both of the CDF and DO Zj cross section measure- 
ments are consistent with the SM expectations at NLO. 
DO searches for anomalous Zjj and ZZj couplings us- 
ing the observed photon Et spectrum in the Zj — >■ ££j 
[Sq l and Zj -)■ i^i^j A LO simulation lia ISTj of 
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Number of events 


W ^ ef 


9.67 ± 


0.30(stat.) 


± 0.48(syst.) 


non- collision 


5.33 ± 


0.39(stat.) 


± 1.91(syst.) 


W/Z + jet 


1.37 ± 


0.26(stat.) 


± 0.91(syst.) 




0.90 ± 


0.07(stat.) 


± 0.12(syst.) 


Total background 


17.3 


± 0.6(stat 


) ± 2.3(syst.) 








33.7 ±3.4 


Nobs 






51 



TABLE IX Summary of background estimates, and the num- 
ber of observed and SM predicted events for the DO analysis 
of the Z'y — J- fvj channel [s^ . 



the Z7 signal is used with NLO corrections \5T) applied 
to the photon Et spectrum. The observed photon Et 
spectrum along with the SM expectation and possible 
anomalous TGC scenarios for two channels is shown in 
Fig. 1161 The Zj — > ££j and Zj — ^ combined lim- 
its on the CP-conserving Z77 and ZZj couplings are 
shown in Fig. [T71 The ID combined limits at 95% CL 



are \hJo\ < 0.033, \hJa\ < 0.0017 
|/ifo| < 0.0017 for A = 1.5 TeV. 



\hlo\ 



< 0.033, and 



4. WW 

Production of W boson pairs in hadron (and lepton) 
collisions involves both the WW'y and WWZ couplings. 
First evidence for W boson pair production was reported 
by CDF using Tevatron Run 1 data [5^ . This process was 
later measured with greater significance by DO and CDF 
using J Cdt = 224-252 pb'^ and J £ dt = 184 pb'^ re- 
spectively, from Run II of the Tevatron [13; [Ml. At LEP, 
WW production has been extensively studied and strin- 
gent limits on anomalous TGC were determined. At the 
Tevatron, much higher WW invariant masses are probed 
compared to LEP because of the higher accessible ener- 
gies. Also, the WW final state is a promising discovery 
channel for the Higgs boson at both the Tevatron and 
the LHC. In hadron collisions, the production of W bo- 
son pairs is most easily observed in the fully leptonic 
decay mode WW — li/ii/. The experimental signature 
of the WW signal in leptonic decay is two isolated high 
Et charged leptons with opposite charge and large 
from the neutrinos. 

Both CDF and DO have measured the WW production 
cross section in fully-leptonic decay and use their data to 
search for anomalous WWj and WWZ couplings. In 
both analyses, the dominant backgrounds are from tt, 
Drell-Yan, other diboson decays, and T^-l-jets where the 
jet fakes an isolated lepton. In the CDF analysis [62l|, ti is 
suppressed by requiring no reconstructed jets with Et > 
15 GeV and \t]\ < 2.5. In the DO analysis the pt 
of the WW system, estimated from the observed charged 
lepton momenta and the ]^ti is required to be small (<20 
GeV/c (ee), 25 GeV/c (e^), or 16 GeV/c(^^)) in order 
to suppress tt decays. 
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FIG. 16 The observed photon Et spectrum along with the 
SM expectation and possible anomalous TGC scenarios for 
(a) Uj ^ and (b) Z-y wy ^ channels. 



The strategy for measuring the cross section for pp — >■ 
WW + X differ between the CDF and DO analyses. In 
the DO analysis, the WW signal yield is determined from 
counting the number of events in excess of the expected 



SM backgrounds using J C dt 
Tab. [X] DO measures 



1.1 fb" 



as shown in 



cr(pp ^ WW + X) 
= 11.5 ± 2.1(stat. 



syst.) ± 0.7(lum.) pb 



[63|iin agreement with the NLO expectation of 12.0 ± 0.7 
pb M- 

In the CDF analysis, the WW signal yield is ex- 
tracted from a fit to the distribution of a matrix element 
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FIG. 17 Two-dimensional bounds (ellipses) at 95% CL on 
CP-conserving (a) ^77 and (b) ZZf couplings. The crosses 
represent the one-dimensional bounds at the 95% CL setting 
all other couplings to zero. The dashed lines indicate the 
unitarity limits for A = 1.5 TeV [H]. 



Process 


Events 


Z/j* (Drell-Yan) 


79.8 ± 18.4 


WZ 


13.8 ± 1.9 


W-y 


91.7 ± 24.8 


W + 1-jet 


112.7 ± 31.2 


ZZ 


20.7 ± 2.8 


ti 


1.3 ± 0.2 


Total Background 


320.0 ± 46.8 




317.6 ± 43.8 


Total Expected 


637.6 ± 73.0 


Data 


654 



TABLE XI Expected number of signal (WW) and back- 
ground events along with the total number of expected and 
observed events in the data for the CDF WW cross section 
measurement 1631. 



Process 


ee 




MM 


Z/'y* — > ee/fifi 


0.27 ± 0.20 


2.52 ± 0.56 


0.76 ±0.36 


Z/j* — > TT 


0.26 ±0.05 


3.67 ±0.46 




ti 


1.10± 0.10 


3.79 ±0.17 


0.22 ±0.04 


WZ 


1.42 ± 0.14 


1.29 ± 0.14 


0.97±0.11 


ZZ 


1.70 ±0.04 


0.09 ±0.01 


0.84 ±0.03 


W-i 


0.23 ±0.16 


5.21 ±2.97 




W-l-jet 


6.09 ± 1.72 


7.50 ± 1.83 


0.12 ± 0.24 


Multijet 


0.01 ±0.01 


0.14± 0.13 




WW 11 


10.98 ±0.59 


39.25 ±0.81 


7.18 ±0.34 


WW It/tt ^ 


W 1.40 ±0.20 


5.18 ±0.29 


0.71 ±0.10 


Total expected 


23.46 ± 1.90 


68.64 ±3.88 


10.79 ±0.58 


Data 


22 


64 


14 



TABLE X Numbers of signal and background events ex- 
pected and number of events observed after the final event 
selection in each channel for the DO WW cross section mea- 
surement 1631. 



likelihood ratio {LRww) discriminant for events using 
J C dt ~ 3.6 fb^^. The events which are fit are required 
to pass relatively loose selection criteria as compared to 
the selection CDF would use for a cross section measure- 
ment based on the event yield alone. Tab. IXII shows the 
expected number of signal and background events along 
with the observed events in the data used to fit for the 
signal. For each event passing the signal selection crite- 
ria, four matrix-element-based (ME) event probabilities 
are calculated corresponding to the production and decay 
processes WW iviv, ZZ Uw, W + 1-jet iv+ 
1-jet, and W7 — £1^ + 7. In the latter two processes, 
the jet or 7 is assumed to have been reconstructed as a 
charged lepton candidate. The event probability for a 
process X is given by 



Px{x) = 



1 



d(j{y) 



e{y) G{x, y) dy 



(7) 



where x represents the observed lepton momenta and 
frp vectors, G{x, y) is a transfer function representing 
the detector resolution, and e{y) is an efficiency function 
parametrized by 77 which quantifies the probability for a 



particle to be reconstructed as a lepton. The differen- 
tial cross section ^^^^ is calculated using leading-order 
matrix elements from the mcfm program |64] and inte- 
grated over all possible true values of the final state par- 
ticle 4-vectors y. The normalization factor (cr) is deter- 
mined from the leading-order cross section and detector 
acceptance for each process. These event probabilities 
are combined into a likelihood ratio 



LRww — 



WW 



'ww + J2i ^jPj 



(8) 



where j = {ZZ, W + 1— jet, W^} and kj is the relative 
fraction of the expected number of events for the j-th pro- 
cess such that kj = 1. The templates of the LRww 
distribution are created for signal and each background 
process given in Tab. IXII 

A binned maximum likelihood is used to extract the 
WW production cross section from the shape and nor- 
malization of the LRww templates. The likelihood is 
formed from the Poisson probabilities of observing Ui 
events in the i-th bin when /i^ are expected. Variations 
corresponding to the systematic uncertainties described 
previously are included as normalization parameters for 
signal and background, constrained by Gaussian terms. 
The likelihood is given by 



n 



Mi 



where 



E 



oik 



1 + /fc^c 



(9) 



(10) 



is the fractional uncertainty for the process k due to 
the systematic c, and Sc is a floating parameter asso- 
ciated with the systematic uncertainty c. The correla- 
tions of systematic uncertainties between processes are 
accounted for in the definition of /^i. The expected num- 
ber of events from process k in the i-th bin is given by 
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FIG. 18 The LRww distributions for the signal [W^W) 
and background processes after a maximum likehhood fit to 
the data for the CDF WW cross section measurement 1621. 



{N^^^)i. The parameter au is an overall normalization 
parameter for process k and is fixed to unity for all pro- 
cesses other than WW, for which it is freely floating. 
The likelihood is maximized with respect to the system- 
atic parameters Sc and aww- The WW cross section 
is then given by the fitted value of aww multiplied by 
a^^'^ipp^WW). 

The fit to the data of the signal and sum of the indi- 
vidually fitted background templates is show in Fig. [TSl 
The measured WW production cross section is 

a{pp WW + X) 

= 12.1 ± 0.9(stat.) L-i.4(syst.) pb 

[6^ in agreement with the NLO theoretical expectation 
of 11.7 ± 0.7 pb m. 

Both of the CDF and DO WW cross section measure- 
ments are consistent with the SM expectations at NLO. 
CDF searches for anomalous WWj and WW Z couplings 
using the observed leading charged lepton pT spectrum 
(see Fig. [in]). In the DO analysis, the subleading (trail- 
ing) lepton pt is also included in a 2D histogram with 
the leading lepton p^ to constrain possible anomalous 
couphngs (see Fig. [20] for the ID projections). 

There are several ways to relate the WWj and WW Z 
couplings in the presence of new physics. This is a con- 
venient prescription to reduce the number of parameters 
since WW production involves both WWj and WW Z 
couphngs. Enforcing SU{2)l <8) C/(l)y symmetry intro- 
duces two relationships between the remaining parame- 
ters: Kz = gi — {kj — l)tan^6'vi/ and Xz — reducing 
the number of free parameters to three 0; Q. Alterna- 
tively, enforcing equality between the WWj and WWZ 



vertices {WWj=WWZ) such that k-, 



and 
two. 



Kz, 



A; 



1 reduces the number of free parameters to 



In the DO analysis, the one-dimensional 95% CL limits 
for A = 2 TeV are determined to be -0.54 < Ak^ < 0.83, 
-0.14 < A^, = Az < 0.18, and -0.14 < A.gf < 



SMW*W" X^ = 0.16 

Background Ag^ = 0.34 

Data Aic, = 0.72 
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FIG. 19 Leading-lepton pr distribution for data compared 
to the SM expectation for the CDF WW — > £iy£u analysis 
[6^ . Also shown is how the expectation would be modified 
by anomalous couplings near the observed limits. 
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FIG. 20 Distributions of the (a) leading and (b) trailing lep- 
ton pt for the DO WW — >■ iuiv analysis [6^. Data are com- 
pared to estimated signal, a(M^M^) — 12 pb, and background 
sum. 



0.30 under the SU{2)l (8) f(l)y-conserving constraints, 



and -0.12 < Ak^ 



Akz < 0.35, with the same A 



limits as above, under the WWj=WWZ constraints. 
One- and two-dimensional 95% CL limits are shown in 
Fig. [5TJ In the CDF analysis, only ID limits on the 
anomalous coupling parameters under the assumption of 
SU{2)l ® U{1)y invariance are reported. The expected 
and observed 95% confidence limits are shown in Tab. lXIll 
where it is evident that the limits are weaker than ex- 
pected. The probability of observing these limits in the 
presence of only SM WW production ranges from 7.1% 
to 7.6% depending on the coupling parameters (Az, Ak^, 
Agf ) and are consistent with a statistical fluctuation of 
SM physics. 



5. WZ 

The W Z final state is not available in e+e~ collisions at 
LEP but can be produced in pp collisions at the Tevatron. 
The study of associated production of a and Z boson 
is important for a number of reasons. The production of 
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FIG. 21 One and two-dimensional 95% CL limits for the DO 
WW ivlv analysis [H when enforcing SU{2)l ® U{1)y 
symmetry at A = 2 TeV, for (a) Ak-^ vs. A^, (b) Aft-y vs. 
Apf , and (c) A^. vs. Agf , each when the third free coupling is 
set to its SM value; limits when enforcing the WW'y=WW Z 
constraints are shown in (d). The curve represents the two- 
dimensional 95% CL contour and the ticks along the axes 
represent the one-dimensional 95% CL limits. An asterisk 
(^) marks the point with the highest likelihood in the two- 
dimensional plane. 





A (TeV) 


A-, — Xz 




Afff 


Expected 
Observed 


1.5 
1.5 


(-0.05,0.07) 
(-0.16,0.16) 


(-0.23,0.31) 
(-0.63,0.72) 


(-0.09,0.17) 
(-0.24,0.34) 


Expected 
Observed 


2.0 
2.0 


(-0.05,0.06) 
(-0.14,0.15) 


(-0.20,0.27) 
(-0.57,0.65) 


(-0.08,0.15) 
(-0.22,0.30) 



TABLE XII Expected and observed limits from the CDF 
WW — >■ li/Eu analysis searching for anomalous TGCs assum- 
ing two different values of the form factor scale A [6^. For 
each coupling limit set, the two other couplings are fixed at 
their SM values. Values of the couplings outside of the given 
observed range are excluded at the 95% confidence level (CL). 



WZ involves the WWZ TGC as shown in the s-channel 
diagram in Fig. [221 Unlike WW production which in- 
volves a combination of the WW^ and WW Z couplings 
such that assumptions regarding their relation must be 
invoked to interpret any anomalies observed in the data, 
the W Z production characteristics can be used to make 
a model-independent test of the SM WWZ coupling. 
Stated slightly differently, in WZ production measure- 
ments, a direct measure of the W Z coupling independent 
of the W^ coupling can be made and compared to the 
SM predictions. The fully-leptonic decay mode of WZ 
provides a clean SM trilepton signal which is analogous 
to the so-called golden mode for discovering supersym- 
metry (SUSY) at the Tevatron via chargino-neutralino 




(b) 



FIG. 22 Leading order (a) t-channel and (b) s-channel dia- 
grams for WZ production at the Tevatron 



production (xf X2) ^^'^ decay. Therefore, an observation 
of the SM WZ trilepton signal represents an important 
experimental milestone in demonstrating sensitivity to 
the SUSY golden mode and other new physics signatures 
in multileptons. 

Prior to the start of Run II at the Tevatron, WZ pro- 
duction had not been observed. The next-to-leading or- 
der (NLO) WZ cross section prediction for pp collisions 
at = 1.96 TeV is 3.7 ± 0.3 pb In October 2006, 
WZ production was first observed by the CDF Collabo- 
ration in the three charged lepton -I- final state using 
1.1 fb~^ of integrated luminosity f65j. The most sensitive 
previous search for WZ production was reported by the 
-DO Collaboration using J C dt = 0.3 fb~^, where three 
WZ — >■ i'viiU candidates were found 66]. The observed 
candidates had a probability of 3.5% to be due to back- 
ground fluctuations, corresponding to a{WZ) < 13.3 pb 
at 95% CL. The DO collaboration published an update to 
[66l | with additional data to measure the WZ production 
cross section and search for anomalous WWZ couplings 

As with other diboson processes, the signal for WZ 
production is most easily measured in fully-leptonic de- 
cay. The experimental signature of WZ production is 
three isolated high Et charged leptons, at least two of 
which having the same flavor and opposite charge with 
invariant mass consistent with decay of a Z boson, and 
large consistent with a neutrino from W decay. To 
observe WZ in fully-leptonic decay, high acceptance for 
charged leptons is required since all three must be de- 
tected to suppress backgrounds from larger cross section 
processes. In the CDF analysis [6^, a novel lepton identi- 
fication strategy was developed to maximize charged lep- 
ton acceptance while keeping the backgrounds compara- 
tively low by exploiting the charge and flavor correlation 
of identified leptons in the events. The standard elec- 
tron and muon identification was combined with forward 
electron candidates beyond the tracking acceptance and 
a "track-only" lepton category consisting of high-quality 
tracks that neither project to the fiducial regions of the 
calorimeters nor are identified as muons by the muon 
chambers. Due to the lack of calorimeter information, 
electrons and muons cannot be reliably differentiated for 
this category, and are therefore treated as having either 
fiavor in the WZ candidate selection. For forward elec- 
trons without a matched track, both charge hypotheses 
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Source 


Expectation ± Stat ± Syst ± Lumi 


Z+jets 


1.21 ± 0.27 ± 0.28 ± - 


ZZ 


0.88 ± 0.01 ± 0.09 ± 0.05 


Zl 


0.44 ± 0.05 ± 0.15 ± 0.03 


tt 


0.12 ± 0.01 ± 0.02 ± 0.01 


Total Background 


2.65 ± 0.28 ± 0.33 ± 0.09 


WZ 


9.75 ± 0.03 ± 0.31 ± 0.59 


Total Expected 


12.41 ± 0.28 ± 0.45 ± 0.67 


Observed 


16 



TABLE XIII Expected number of events in the signal re- 
gion for W Z and the background contributions for the CDF 
WZ — >■ Illy analysis [g^. "Lumi" refers to the integrated lu- 
minosity uncertainty, which is absent for the Z-|-jets because 
it is determined from the same data set. 



are considered when forming W Z candidates, since the 
charge is determined from the track curvature. 

Other SM processes that can lead to three high-p^ lep- 
tons include dileptons from the Drell-Yan Z/^* process 
(DY), with an additional lepton from a photon conver- 
sion {Z^) or a misidentified jet (Z-|-jets) in the event; 
ZZ production where only three leptons are identified 
and the unobserved lepton results in $rp-^ and a small 
contribution from ti — >■ WbWb, where two charged lep- 
tons result from the W boson decays and one or more 
from decay of the 6-quarks. Except for tt, these back- 
grounds are suppressed by requiring > 25 GeV in the 
event, consistent with the unobserved neutrino from the 
leptonic decay oi aW boson. 

Using J C dt = l.l fb~^, 16 events with an expected 
background of 2.7 ± 0.4 events are observed in the CDF 
analysis, as shown in Fig. IXIIII This corresponds to a 
6.0(7 observation of WZ production when $!rp shape in- 
formation is included. The measured cross section is 

a{pp ^ WZ + X)^ 5.0tli (stat. + syst.) pb, 

consistent with the NLO expectation. Fig.l^shows some 
important kinematic distributions for the 16 WZ candi- 
dates, which are in good agreement with the SM. 

In the DO analysis [67|, events with three recon- 
structed, isolated charged leptons (electrons or muons) 
and J^rj. > 20 GeV were used to measure the WZ produc- 
tion cross section and search for anomalous WW Z cou- 
plings using the observed Z boson (Pt) distribution. 
Each W Z candidate event must contain a like-flavor lep- 
ton pair with invariant mass close to the Z boson mass. 
For eee and decay channels, the lepton pair with 
invariant mass closest to that of the Z boson mass is 
chosen to define the Z boson daughter particles. Using 
/£ = 1.0 fb"\ a total of 13 WZ candidates are ob- 
served with 4.5 ± 0.6 expected background events and 
9.2 ± 1.0 expected WZ signal events, corresponding to 
a 3.0(7 signal significance. The breakdown by trilepton 
flavor classification for both the CDF and DO analyses is 
shown in Tab. IXIVI Fig. [24] shows versus the dilcp- 
ton invariant mass for the background, the expected W Z 



Flavor CDF Analysis DO Analysis 



Classification 


Expected 


Data Expected Data 


e e 6 


2.7 ± 0.2 


6 3.5 ± 0.2 2 




2.0 ± 0.2 


9 7 + 9 1 


e fj, fj, 


1.5 ± 0.1 


1 4.2 ± 0.5 8 


fi fj, fi 


1.2 ± 0.1 


1 3.4 ± 0.4 2 


e e It 


2.0 ± 0.2 


5 - - 


e fi It 


1.3 ± 0.1 


2 - - 


fJ- fJ. It 


1.1 ± 0.1 


1 - - 


e It It 


0.5 ± 0.1 


- - 


filt it 


0.2 ± 0.1 


- - 


TABLE XIV Summary of the 


expected and observed yields 


for the CDF Ml 


and DO l&E 


WZ lllv analyses. In the 


classification column. It denotes a track-only lepton candidate 


having unknown flavor which 


is only relevant for the CDF 


analysis. 






A = 1.5 TcV 


A = 2.0 TeV 


-0.18 < Xz 


< 0.22 


-0.17 <\z < 0.21 


-0.15 < Agf < 0.35 


-0.14 < Afff < 0.34 


-0.14 < Akz = 


Agf < 0.31 


-0.12 < Akz = Apf < 0.29 



TABLE XV One-dimensional 95% CL intervals on Az, Agf , 
and Akz for two sets of form factor scale. A, for the DO 
WZ mu analysis [gj]. 



signal, and the data, including the candidates for the DO 
analysis. The W Z production cross section measured by 
DO is 

f7(pp WZ + V) = 2.7li;3 (stat. + syst.) pb, 

where the ±1ct uncertainties are the 68% CL limits from 
the minimum of the negative log likelihood. The uncer- 
tainty is dominated by the statistics of the number of 
observed events. 

The three CP-conserving WW Z coupling parameters 
\z, Agf , and Ak^ are constrained in the DO analysis by 
comparing the measured cross section and distribu- 
tion to models with anomalous couplings. A comparison 
of the observed Z boson distribution with predictions 
from Monte Carlo simulation is shown in Fig. [55] Tab. 
[XV] presents the one-dimensional 95% CL limits on A^, 
Agf and Ak^. Fig.[2n]presents the two-dimensional 95% 
CL limits under the assumption Ayf = Ak^ for A = 2 
TeV. 



6. ZZ 

The production of Z pairs is predicted within the SM 
to have the smallest cross section among the diboson pro- 
cesses. It has been been observed in e+e^ collisions at 
LEP [68] , but not in hadron collisions as of the start of the 
Tcvatron Run IL As a window to new physics, ZZ pro- 
duction is particularly interesting because of the absence 
of ZZ^i and ZZZ couphngs in the SM (see Fig.]!?]), and 
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FIG. 23 Distributions for WZ candidates of (a) the (b) the dilepton invariant mass for the same-fiavor opposite-sign 
dilepton pair closest to the Z mass, and (c) the W transverse mass calculated from the remaining lepton and the l^j,, for the 
CDF WZ — >■ iUv analysis In (a) and (b), the arrows indicate the signal region. 
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FIG. 24 versus dilepton invariant mass of W Z ^ Mif 
candidate events in the DO WZ analvsis[67j. The open boxes 
represent the expected WZ signal. The gray boxes represent 
the sum of the estimated backgrounds. The black stars are 
the data that survive all selection criteria. The open circles 
are data that fail either the dilepton invariant mass criterion 
or have $j, < 20 GeV. 
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FIG. 25 The reconstructed Z boson pt of the WZ candi- 
date events for the DO WZ -j- mu analysis The solid 
histogram is the expected sum of signal and background for 
the case of the WW Z coupling parameters set to their SM 
values. The dotted and double dotted histograms are the ex- 
pected sums of signal and background for two different cases 
of anomalous WW Z coupling parameter values. 



because of the very low backgrounds in the four charged- 
lepton channel Higgs boson decay can contribute to ZZ 
production; hovifever, this channel is generally not com- 
petitive With H VFVF^*^ as a discovery channel at the 
Tevatron collision energy and integrated luminosity. 

As is the case in for WW and WZ production, the 
ZZ state is most easily observed in the fully leptonic 
mode at a hadron collider. The ZZ — > £U£ process 
is rare but predicted to be nearly background free in 
the SM, with Z-l-jets (jets reconstructed as charged lep- 
tons) as the only non-negligible background. Having 
large total charged lepton acceptance in the experiments 



is crucial due to the high lepton multiplicity in the fi- 
nal state. The ZZ — process can also contribute 
sensitivity to the search for ZZ production although it 
suffers from large continuum WW backgrounds. The 
fuU SM process is pp Z/j* Z/'y*, where the two 
Z/7* interfere with one another. For brevity, we denote 
Z/7* Z/j* as ZZ throughout and indicate the dilepton 
invariant mass rangc(s), where applicable. The next-to- 
leading order (NLO) ZZ cross section for pp collisions at 
^/s ~ 1.96 TeV is 1.4 ± 0.1 pb in the zero-width Z boson 
approximation [64l |. 

Using J C dt = 1 fb~^, the DO Collaboration searched 
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FIG. 26 Two-dimensional 95% CL contour limit in Agf — 
Akz versus AAz space (inner contour) for the DO WZ — > 
lUv analysis [bJ]- The form factor scale for this contour is 
A = 2 TeV. The physically allowed region (unitarity limit) is 
bounded by the outer contour. The cross hairs are the 95% 
CL one-dimensional limits. 
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FIG. 28 Limits on anomalous couplings for A = 1.2 TeV for 
the DO ZZ MM analysis (a) fJo vs. /|,, (b) f2o vs. 

/sO' (c) fio vs. /so, and (d) f^g vs. fio, assuming in each 
case that the other two couplings are zero. The inner and 
outer curves are the 95% CL two-degree of freedom exclusion 
contour and the constraint from the unitarity condition, re- 
spectively. The inner cross hairs are the 95% CL one-degree 
of freedom exclusion limits. 




FIG. 27 Leading order diagram for ZZ production at the 
Tevatron. The s-channel diagram involving the neutral TGC 
is absent in the SM. 



for ZZ — T' ££i£ production and set a limit on the cross 
section of a{ZZ) < 4.4 pb at 95% CL ^]. They con- 
strain possible ZZj and ZZZ couplings based on the 
observed four lepton yield with the added requirement 
that AI{££) > 50 (70) for electrons (muons). The value of 
the cut was chosen based on the dilepton invariant mass 
resolution. This lower M{££) requirement was included 
because the Monte Carlo generator [i^l used to constrain 
possible anomalous couplings does not include contribu- 
tions from off-shell Z bosons. Without the M{££) cut, 
one /i/iee event is observed in the data. This event is 
removed by the M{££) cut used to constrain anomalous 
couplings. ID and 2D limits on anomalous ZZj and 
ZZZ couplings are determined using A = 1.2 TeV. The 
95% CL ID limits are -0.28 < fg < 0.28, -0.26 < fJo < 
0.26, -0.31 < fg < 0.29, and -0.30 < /g^ < 0.28. The 
95% CL 2D contours fJo vs. fg, f2o vs. /50, /fo vs. fg, 
and /gg vs. are shown in Fig. [28l 

The CDF Collaboration uses /£ = 1.9 fb"^ to 
search for ZZ production in a combination of the ZZ — >■ 
££££ and ZZ ££vv channels 7^. To maximize the 



acceptance, lepton candidates are constructed out of all 
reconstructed tracks and energy clusters in the EM sec- 
tion of the calorimeter. This is done with the same lepton 
identification criteria used in a previous CDF measure- 
ment of WZ production The ZZ ££££ candidates 
are selected from events with exactly four charged-lepton 
candidates and at least two same-flavor, opposite-sign 
lepton pairs are required for the event to be accepted. 
As in the W Z analysis, trackless electrons are consid- 
ered to have either charge, and track-only leptons either 
flavor. One pair must have invariant mass M^+i- in the 
range [76, 106] GeV/c^, while the requirement for the 
other pair is extended to [40, 140] GeN/c to mcrease 
the acceptance for off-shell Z decays. 

The ZZ — >■ ££££ candidate events are separated into 
two exclusive categories based on whether or not they 
contain at least one forward electron without a track. 
This is done because the background from Z7+jets is 
much larger in candidates with a forward trackless elec- 
tron. The expected signal, expected background, and 
observed yields are shown in Tab. IXVII 

The ZZ — ££vv candidates are selected from events 
with exactly two oppositely-charged lepton candidates 
excluding events with forward electrons without a track 
which are contaminated by large VF7 backgrounds. Aside 
from ZZ production, other SM processes that can lead 
to two high-pT leptons include events from DY, a W 
decay with photon (VF7) or jet (W^-|-jets) misidentified 



23 





Candidates without a 


Candidates with a 


Category 


trackless electron 


trackless electron 


zz 


1.990 ± 0.013 ± 0.210 


0.278 ± 0.005 ± 0.029 


Z+jets 


0.014t;3 [;j? ± 0.003 


0.082t;3 !5|g ± 0.016 


Total 


2.004to.°i^ ± 0.210 


0.3601;^ E;^^ ± 0.033 


Observed 


2 


1 



TABLE XVI Expected and observed number of ZZ UU 
candidate events for the CDF ZZ analysis [t^]. The first 
uncertainty is statistical and the second one is systematic. 
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FIG. 29 Distribution of the discriminating variable logj^Q(l 
LR) for the CDF ZZ Uvv search 



as a lepton; and tt, WW, and WZ production. 

There are 276 events after the event selection (which 
contains a speciahzed high l^rp to suppress primarily 
DY) of which only 14 ± 2 are expected to be from the 
ZZ tlvv process in the SM. Approximately half of the 
yield is expected to be due to the WW process. However, 
ZZ — )■ Uvv and WW have different kinematic properties 
which are exploited to statistically separate the contri- 
bution of these two processes to the data. The approach 
used by CDF is identical to that used in the WW cross 
section measurement described in Section IIII.E.4I An 
event- by- event probability density is calculated for the 
observed lepton momenta and $rp using leading order cal- 
culations of the differential decay rate for the processes 
[63 | . A likelihood ratio discriminant LR is formed which 
is the signal probability divided by the sum of signal and 
background probabilities LR — Pzz /{Pzz + Pww ) ■ The 
distribution of logiQ(l — LR) for the data compared to 
the summed signal and background expectation is shown 
in Fig. [m 
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FIG. 30 Distribution of four lepton invariant mass in data, 
expected signal, and expected background for the DO ZZ — >■ 
UU analysis 



The ZZ UU and ZZ — ^ Uvv results were combined 
through a likelihood fit which includes the two 4-lepton 
bins and the log]^o(l ^ distribution for the dilep- 
tons. The p- value for the ZZ — )■ Uvv alone is 0.12 and 
the combined p-value is 5.1 x 10~^ corresponding to a 
significance equivalent to 4.4 standard deviations. The 
ZZ cross-section is obtained by fitting the data for the 
fraction of the expected SM yield in the full acceptance 
and scaling the zero- width Z boson approximation cross- 
section by that fraction. The measured cross section is 
a{pp ZZ) = 1.4l° g pb, consistent with the SM ex- 
pectation. This is the first evidence of a ZZ signal with 
greater than 4cr significance in hadron collisions. 

The DO Collaboration uses J C dt — l.Tfb^"'^ to search 
for ZZ production in the ZZ UU channel This 
analysis follows up the analysis from [69l | using more data 
and tighter dilepton invariant mass requirements but 
does not explicitly search for anomalous couplings. Four 
lepton events are selected using identified muons with 
|?7| < 2 and electrons that are either central (jry < 1.1) 
or forward 1.5 < |?7| < 3.2. It is required that one pair 
have the same flavor with invariant mass > 70 GeV/c 
and another pair have invariant mass > 50 GeV/c^. 

Tab. IXVIII summarizes the expected signal and back- 
ground contributions to each subchannel, as well as the 
number of candidate events in data. The total signal 
and background expectations are 1.89 ± 0.08 events and 
0.14tn n9 events, respectively. A total of three candidate 



-0.02 

events is observed, with two in the Ae^c subchannel ( "C" 
refers to the number of central leptons) and one in the 
4/i subchannel. Fig. [501 shows the distribution of the four 
lepton invariant mass for data and for the expected signal 
and background. 

Using a log-likelihood ratio test statistic of the yields 
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Subchannel 


4e2c 


4e3c 


4:64,0 


4fi 


2/i2eoc 


2^2610 


2^2620 


Luminosity (fb^'^) 


1.75 ±0.11 


1.75 ±0.11 


1.75 ±0.11 


1.68 ±0.10 


1.68 ±0.10 


1.68 ±0.10 


1.68 ±0.10 


Signal 


0.084 ± 0.008 


0.173 ±0.015 


0.140 ±0.012 


0.534 ± 0.043 


o^s+''-°°'^ 

U.UOO_Q QQg 


0.352 ± 0.040 


U.OOO_Q Q44 


Z(7)+jets 


U.UOU_Q QQg 


o.oi8«:r 




0.0003 ± 0.0001 


n3+"°2 


0.05 ±0.01 


0.008_o („)3 


tt 










0012+°"''^'^ 


0.005 ± 0.002 


0007+" °°°^ 

U.UUUr^Q QQQ5 


Observed events 








2 


1 












TABLE XVII The integrated luminosity, expected number of signal {Z/"/* Z/j*) and background events (tt and Z(7)±jets 
which includes aU Vy/^/7+jets contributions), and the number of observed candidates in the seven ZZ — ^ £^£~£ '^l ~ sub- 
channels for the DO ZZ — >■ Mi analysis [tJ] . Uncertainties reflect statistical and systematic contributions added in quadrature. 



and pseudo-experiments, the p-value is determined to be 
4.3xl0~^ which corresponds to a 5.3 a observed signif- 
icance (3.7(7 expected). The measured cross section in 
the ZZ -J> mi channel is 

a{ZZ + X) = l.l^tl^ (stat.) ± 0.13 (syst.) pb. 

This result is combined with the results from an in- 
dependent ZZ Uvv search [tI], and the previous 
ZZ — )• Ull analysis {i^l taking into account the system- 
atic uncertainty correlations between subchannels and 
among analyses. The resulting p- value is 6.2xl0~^, and 
the significance for observation of ZZ production in- 
creases to h.l<7 (4.8cr expected). This is the first ob- 
servation of ZZ production at a hadron collider. The 
combined cross section is 

a{ZZ ± X) = 1.60 ± 0.63 (stat.)+";}7 (syst.) pb, 

consistent with the SM expectation. 

7. WV {y = W, Z) 

The production of vector boson pairs {WW, WZ, 
and ZZ) have been observed at the Tevatron in decay 
modes where both vector bosons decay leptonically. The 
semileptonic decay modes where one of the vector bo- 
son decay hadronically has larger branching fraction as 
compared to the fully leptonic modes but significantly 
larger backgrounds from jets produced in association 
with a W boson. As a result, simple event counting 
above background cannot be used to observe dibosons in 
semileptonic decay at the Tevatron and advanced anal- 
ysis techniques utilizing multivariate event classification 
are required to statistically separate signal from back- 
ground. Analysis of the £vjj final state provides an ex- 
cellent testbed for such advanced techniques to extract 
small signals from large backgrounds in real hadron col- 
lision data that has great relevance to Higgs boson (e.g. 
WH — ^ ivbb) and new physics searches in final states 
involving jets. In addition, semileptonic diboson decay 
can provide a more sensitive search for anomalous tri- 
linear gauge couplings than fully leptonic decay modes 
since anomalous TGCs enhance production at high gauge 



boson momentum where the signal-to-background ratio 
improves. 

The dijet mass resolution of the CDF and DO detectors 
is not good enough to distinguish hadronically decaying 
W bosons and Z bosons. As a result, the search for 
diboson production in the ii'jj final state is a search for 
the sum of WW and WZ production. The experimental 
signature of the WW + WZ signal in semileptonic decay 
is one isolated high Et charged lepton, large from 
the neutrino produced in W decay, and at least two high 
Et jets. 

In principle, one could search for SM WZ + ZZ pro- 
duction in the jj££ channel. Large Z+jets backgrounds 
makes this a very difficult channel to observe the SM sig- 
nal, although the use of 6-tagging can improve the sensi- 
tivity. The use of the bb£i channel to search for associated 
Higgs boson production is described in Section IV. B. 21 

Both DO and CDF have searched for WW + WZ pro- 
duction and anomalous WWj, WWZ couplings in the 
ivjj final state (73l - [76l |. In these analyses, Il^+jets is the 
dominant background. Other significant backgrounds in- 
clude Z+jets, tt, single top quark, and QCD multijet 
production. Backgrounds are suppressed by requiring 
a minimum W transverse mass Mt{W) in each event, 
which is the transverse mass of the charged lepton and 
system. 

The strategy for extracting the semileptonic WW + 
WZ signal yield after basic event selection differs be- 
tween the DO and CDF analyses. In the DO analysis 
thirteen kinematic variables (e.g. dijet mass) demon- 
strating a sensitivity to distinguish signal and back- 
ground are used as input to a Random Forest (RF) [77| 
multivariate event classifier using J C dt = 1.1 fb""'^. 
Fig. [31] shows the RF output distribution after a fit to 
the WV signal and W+jets background contributions. 
Fig. [22] shows the dijet mass distribution using the re- 
sults of the RF output fit. The dominant systematic un- 
certainties arise from the modeling of the M^+jets back- 
ground and the jet energy scale. The probability for the 
background to fiuctuate to give an excess as large as that 
observed in the data is < 5.4 x 10~^, corresponding to a 
4.4(7 signal significance. This result is the first evidence 
for WV production is lepton + jets events at a hadron 
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FIG. 31 For the DO analysis (a) the RF output distri- 
bution from the combined evqq and ^ii^qq channels for data 
and simulated event predictions following the fit of simulation 
to data, (b) A comparison of the extracted signal (filled his- 
togram) to background-subtracted data (points), along with 
the ±1 standard deviation (s.d.) systematic uncertainty on 
the background. The residual distance between the data 
points and the extracted signal, divided by the total uncer- 
tainty, is given at the bottom. 



FIG. 32 For the DO analysis [t^I, (a) the dijet mass distri- 
bution from the combined evqq and fivqq channels for data 
and simulated event predictions following the fit to the RF 
output, (b) A comparison of the extracted signal (filled his- 
togram) to background-subtracted data (points), along with 
the ±1 standard deviation (s.d.) systematic uncertainty on 
the background. The residual distance between the data points 
and the extracted signal, divided by the total uncertainty, is 
given at the bottom. 



collider. The measured cross section is 

a{WV + X) = 20.2 ± 4.5(stat. + syst.) pb 

consistent with the NLO SM cross section of 16.1 ± 0.9 
pb [6|. 

In the CDF analysis [T^I, two different methods are 
used to extract the WW + WZ signal from the data. 
In the first method ("dijet method"), the dijet invariant 
mass Mjj is used to extract a signal peak from data cor- 
responding to J C dt = 3.9 fb~^. The second method 
("ME method") takes advantage of more kinematic in- 
formation in the event by constructing a discriminant 
based on calculations of the differential cross sections of 
the signal and background processes data corresponding 
to J C dt — 2.7 fb~^. The ME method was discussed in 
Section IIILE.4I where it was used in the WW — ?> iviv 
analysis and has been used in a number of different CDF 
and DO analyses. 

The normalization of the Z+jets background is based 
on the measured cross section [78| . The tt and single top 
background normalizations are from the NLO predicted 



cross sections [t^; H^. The efficiencies for the Z-|-jets, 
tt, and single top backgrounds are estimated from sim- 
ulation. The normalization of the multijet background 
is estimated by fitting the J^j, spectrum in data to the 
sum of all contributing processes, where the multijet and 
W^-|-jets normalizations float in the fit. In the final signal 
extractions from both methods, the multijet background 
is Gaussian constrained to the result of this fit and 
the W^-f jets background is left unconstrained. 

In the dijet method, the signal fraction in the data is 
estimated by performing a fit to the dijet invariant 
mass spectrum, separately for electron and muon events. 
Fig. l33l shows the dijet invariant mass distribution of data 
compared to the fitted signal and background contribu- 
tions. The observed (expected) signal significance using 
the dijet method is 4.6ct (4.9cr) for electrons and muons 
combined. The measured cross section is 

a{WV + X) = 14.4 ± 3.1(stat.) ± 2.2(syst.) pb 

In the ME method, the calculated event probabili- 
ties are combined into an event probability discriminant 
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FIG. 33 For the CDF analysis [7J|, the (a) dijet invariant 
mass distribution of reconstructed W/Z jj candidates 
compared to the fitted signal and background components, 
and (b) corresponding background subtracted distribution. 



EPD = Ps^gnal/{Ps^g nal T t 



background 



), where P. 



signal 



Pww + Pwz and Pha ckground — Pw+jets ~t~ Psingle top- 

Templates of tfie EPD generated for all signal and back- 
ground processes are used in a binned likelihood fit for 
tlie signal yield observed in the data, as shown in Fig.lMl 
Fig. inS] shows the dijet mass in bins of EPD, where 
it can be seen that the low EPD bin contains very 
little signal as compared to background. Events with 
EPD > 0.25 have a dijet mass peak close to the expected 
W/Z mass, and the signal-to-background ratio improves 
with increasing EPD. The observed (expected) signal 
significance using the ME method is 5.4tT (5. la) which 
represents the first observation of WW + WZ production 
in the lepton -I- jets channel. The measured cross section 
using the ME method is 

aiWV + X)^ 17.7 ± 3.1(stat.) ± 2.4(syst.) pb. 

In the Mjj method the largest systematic uncertainties 
are due to the modeling of the EW and multijet shapes, 
about 8% and 6% respectively. In the ME method the un- 
certainty in the jet energy scale is the largest systematic 
uncertainty, at about 10%, which includes contributions 
both from the signal acceptance and from the shapes of 
the signal templates. In the Mjj method this uncertainty 
is about 6%. The combined cross section for the dijet 
and ME methods, with consideration of statistical and 
systematic uncertainties, is 

cr{WV + X) = 16.0 ± 3.3(stat. + syst.) pb. 

Both of the CDF and DO WW+WZ cross section mea- 
surements are consistent with the SM NLO expectations. 
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FIG. 34 Observed EPD distribution superimposed on distri- 
bution expected from simulated processes for the CDF anal- 
ysis [zl] 
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FIG. 35 For the CDF analysis [TJ]. Mjj for events with (a) 
EPD < 0.25, (b) 0.25 < EPD < 0.5, (c) 0.5 < EPD < 0.75, 
and (d) EPD > 0.75. 



In [75(], CDF searches for anomalous WWj and WW Z 
couplings using the observed spectrum of the charged lep- 
ton pt from a W decay and J C dt — 0.35 fb^^. Fig. 
shows good consistency between the SM expectation for 
W boson pt and the data which is used to constrain 
possible anomalous couplings. In order to increase the 
sensitivity to anomalous couplings, these data are com- 
bined with the photon Et spectrum from the Wj anal- 
ysis [4^ described in Section UlI.E. 21 that constrains pos- 
sible anomalous WWj couplings. Tab. IXVIIII summa- 
rizes the resulting anomalous coupling limits under the 
assumption of equal WWj and WW Z couplings. 

The DO search for anomalous couplings in the Ivj] 
channel [tI] is based on the same data that was used 
to obtain the first evidence for semileptonic decays of 
WW/WZ boson pairs in hadron collisions [t^. In con- 
trast to the CDF analysis just described, the DO analysis 
uses the dijet px information to place limits on WWj/Z 
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FIG. 36 The pj!' distribution used for setting the anomalous 
TGC limits in the CDF analysis (z!]. Shown are the SM ex- 
pectation, data points, and one anomalous coupling scenario 
with A = Ak = 0.5 and A = 1.5 TeV. 





A 


A/t 




(-0.28, 0.28) 


(-0.50, 0.43) 




(-0.21, 0.19) 


(-0.74, 0.73) 


Combined 


(-0.18, 0.17) 


(-0.46, 0.39) 



TABLE XVIII Allowed anomalous coupling ranges for A — 
1.5 TeV at 95% CL, fixing the other coupling to the SM value 
and assuming equal WWy and WW Z couplings for the CDF 
analysis [tsI ]. 



anomalous couplings. Fig. l37l shows the dijet pt distribu- 
tion after a fit for the individual signal and background 
contributions. No excess at high dijet pT that would 
suggest that non-standard couplings are present is evi- 
dent in the data. The one-dimensional 95% CL limits for 
the WW^/Z anomalous coupling parameters under sepa- 
rate assumptions of SU{2)l®U{1)y symmetry and equal 
WW-i and WWZ couphngs are shown in Tab.|XlXl The 
2D limits under the assumption of SU (2)l ® U{1)y sym- 
metry [63 is shown in Fig. |38] and for equal coupling as- 
sumption they are shown in Fig. 1391 



8. VV (V = W,Z) 

In Sections lIII.E.4llnLE.6[ we summarized the observa- 
tions of heavy vector boson pairs {WW, WZ, ZZ) when 
each boson decays to one or more charged leptons. In 
Section IIII.E.7[ the measurements of a hadronically de- 
caying heavy vector boson produced in association with 
a W boson that decays to a charged lepton and a neu- 
trino was described. Observing diboson production in 
semileptonic decay mode required sophisticated analy- 
sis techniques because of overwhelming M^-fjets in addi- 
tion to the presence of a WW + WZ signal. The ap- 
proach of searching for diboson production in an more 
inclusive way can be taken a step further by dropping 
the requirement of a reconstructed charged lepton and 



Data 
■ Diboson Signal 

□ W-njets 

□ Z-f-jets 

□ Multijet 




100 150 
(a) 



300 



Dijet (GeV/c) 



D0 Run II, 1 . 1 fb' 4- Data - (SM Sig-hBkgd) 
95% CL for K 
---95%CLforA. 



® 

a 2 



Ir 



g 
PJ 



-Ir 
-2r 



50 



100 



150 
(b) 



200 



250 



300 



Dijet (GeV) 



FIG. 37 For the DO analysis [zil, the (a) dijet pt distribu- 
tion of combined (electron-|-muon) channels for data and SM 
predictions following the fit of simulated events to data and 
(b) difference between data and simulation divided by the 
uncertainty (statistical and systematic) for the dijet pr dis- 
tribution. Also shown are the simulated signals for anomalous 
couplings corresponding to the 95% CL limits for An and A in 
the LEP parametrization scenario. The full error bars on the 
data points reflect the total (statistical and systematic) un- 
certainty, with the ticks indicating the contribution due only 
to the statistical uncertainty. 



instead looking at the $Jrp + jets final state. This ap- 
proach includes an additional contribution from Z ^ vv 
over the WW/WZ — ^ (i/jj mode so that it measures 
WW + WZ + ZZ production and allows for an event 
with a diboson decay to be accepted even if one or more 
charged leptons fail a selection cut. The experimental 
challenge is that this mode requires a tighter $!rp cut as 
compared to WW/WZ — >■ iiyjj and a better overall un- 
derstanding of high reconstructed $!rp tails from low true 
]/]j< events with large cross section like QCD multijets. 

The CDF CoUaboration uses J C dt ^ 3.5 fb"^ to 
search for the production of heavy vector boson pairs 
{VV, V=W, Z) where one boson decays to a dijet final 
state [8l| . The most significant backgrounds to the dibo- 
son signal are W{tv)+]eis, Z{vv)+]eis, and QCD mul- 
tijet production. Other less significant backgrounds in- 
clude Z{U)+iets, tt, and single i-quark production. The 
majority of events are collected using an inclusive 1/Jrp 
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SU{2)l ® U{1)y symmetry 
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^ _ 1 07+0.26 
'^1 — J^-"' -0.29 

1 n/l+0.11 


A = 
A = 


n nn+° °® 





95% CL 






A A-y Az^ 


Agf 


SU{2)l 0U{1)y symmetry 


-0.44 < 


< 0.55 


-0.10 < A < 0.11 


-0.12 < Agf < 0.20 


Equal couplings 


-0.16 < Ak 


< 0.23 


-0.11 < A < 0.11 





TABLE XIX The most probable values with total uncertainties (statistical and systematic) at 68% CL for k^, A, and along 
with observed 95% CL one-parameter limits on Ak-,, A, and Agf measured in 1.1 fb~^ of WW/WZ Ivjj events with Ajvp — 
2 TeV for the DO analysis 



trigger which requires $j. > 45 GeV. In the offline se- 
lection, events are required to have If^rp > 60 GeV and 
exactly two reconstructed jets with Et > 25 GeV and 
\V\ < 2.0. 

The shape and normalization of the multijet back- 
ground are determined from the data. A tracking- 
based missing momentum vector, ^j,, analogous to the 

calorimeter-based l^j,, is constructed from the vector sum 
of the transverse momenta of particles measured in the 
tracking system, and is largely uncorrelated to fj, for 
events where jets are misreconstructed. In the absence 
of $!rp arising from mismeasurement in the calorimeter, 

the ]prp and ^j, will be aligned in most events. At large 

values of A<j>{l^rp,^j,), multijets are expected to be the 
dominant component of the data. The dijet mass shape 
and normalization for the multijet background that re- 
mains after the event selection is determined by select- 
ing events with A(f){]^j,,^rj^)>1.0 and subtracting out the 
non-multijet backgrounds. The normalization is scaled 
up to account for the multijet background contamina- 
tion in the region A0(^y,|(j,)<l.O. The shape of the 
multijet background is fit to an exponential in Mjj to 
derive a dijet mass template for use in the Mjj unbinned 
extended likelihood fit performed to extract the diboson 
signal. The distribution of A^-IS* observed in the data 

is in good agreement with the expectations as shown in 
Fig. 1301 giving confidence in the validity of the multijet 
background model. 

Three Mjj template distributions are used in the fit: 
the first is F-|-jets and t-quark production and is taken 
from Monte Carlo simulation; the second is the multijet 
template, where the slope and normalization are Gaus- 
sian constrained to their previously measured values; the 
third template describes the signal. The signal shape is 
comprised of the WW, WZ, and ZZ distributions ob- 
tained from a Gaussian -t- polynomial fit to the signal 
Monte Carlo simulation where the mean and the width 
of the Gaussian distribution are linearly dependent on the 
jet energy scale (JES). The uncertainty associated with 
the JES is the dominant source of systematic uncertainty 
of the diboson cross section measurement. 

Fig. |4T] shows the fit result and a comparison between 
the expected signal and data after background subtrac- 



tion. The signal significance is reported to be greater 
than 5.3(7 from the background-only hypothesis. This 
represents the first observation in hadronic collisions of 
the production of weak gauge boson pairs where one bo- 
son decays to a dijet final state. The measured cross 
section is 

cr{WW + WZ + ZZ + X) 

= 18.0 ± 2.8(stat.) ± 2.4(syst.) ± l.l(lumi.) pb, 

consistent with the SM NLO prediction of 16.8 ± 0.5 pb 




9. Cross Section and Anomalous TGC Summary 

A summary of the diboson production cross sections 
at the Tevatron with comparisons to the theory expecta- 
tions can be found in Tab. IXXl and Fig. |42] A summary 
of limits on anomalous charged and neutral TGCs at the 
Tevatron can be found in Tab. IXXj and Tab. IXXlTl re- 
spectively. 

The experimental progress on diboson physics during 
Run II at the Fermilab Tevatron has been truly remark- 
able. This includes the first observation of WZ produc- 
tion and the first observation at a hadron collider of ZZ 
production and of weak gauge boson pair production in 
semileptonic decay. The associated diboson cross section 
measurements (many of which are significantly improved 
from Run I) and detailed studies of their kinematic prop- 
erties provide stringent tests standard electroweak theory 
and constraints on possible new physics coupling to the 
EW gauge bosons. 

IV. TOP QUARK 

A. Top Quark Mass 

As discussed in section fV.Ai constraints on the Higgs 
boson mass from fits to data assuming the SM are limited 
by the experimental precision of the W boson and top 
quark masses. Measurements of the top quark mass thus 
have direct implications for the SM, and these measure- 
ments are a central component of the programs for both 
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Process 


Channel (fb"') 




Cross section (pb) 


Expt. 




Theory 


Measured 


Reference 



pp _^ -y + X — >• £i^7 + X 
















(El > 7 GeV ARf ^ > 71 




~ 0.2 


19.3 ± 1.4 


[43] 


18.1 ± 3.1 (stat. ± syst.) 


CDF 


[461 


{E"r[, > 8 GeV, ARe,^ > 0.7) 




~ 0.2 


16.0 ± 0.4 


[431 


14.8 ± 1.9 (stat. + syst.) ± 1.0 (lumi.) 


DO 


[471 


pp _^ + X ^^7 + X 
















{e:^ > 7 GeV, Ai?^,^, > 0.7, 
















Met > 40 GeV/c") 


U-y 


~ 0.2 


4.5 ±0.3 


[57] 


4.6 ± 0.6 (stat. ± syst.) 


CDF 


[46] 


(EJ > 7 GeV, ARi^^ > 0.7, 
















M« > 30 GeV/c^) 




~ 1.1 


4.7 ±0.2 


[57] 


5.0 ± 0.3 (stat. ± syst.) ± 0.3 (lumi.) 


DO 


[56] 


— !> + X ^ i/vy + X 
















> 90 GeV) 


vvy 


3.6 


0.039 ± 0.04 


[57] 


0.032 ± 0.09 (stat. ± syst.) ± 0.02 (lumi.) 


DO 


[5^ 


pp WW + X 




3.6 


11.7 ±0.7 


[64] 


12.9 ±0.9 (stat.) + !-^ (syst.) 

\ / — 1.4 V ■/ / 


CDF 


[621 




Mp 


~ 1.1 






11.5 ± 2.1 (stat. ± syst.) ± 0.7 (lumi.) 

5.0ti:^ (stat.) ± 0.4 (syst.) 
2.7li:3 (stat. ± syst.) 


DO 


[63] 


pp^WZ + X 


mv 

mv 


1.1 
1.0 


3.7 ±0.3 


[64] 


CDF 
DO 


[65] 
[67] 


pp^ ZZ + X 


mi, iivv 


1.9 


1.4 ±0.1 


[64] 


l-4lo;6 (stat. + syst.) 


CDF 


[70] 




im, iivv 


1.7 






1.6 ± 0.6 (stat.) ± 0.2 (syst.) 


DO 


[7JJ 


pp WV + X {V -^W,Z) 


Ivqq 


2.7 


16.1 ±0.9 


[64] 


16.0 ± 3.3(stat. ± syst.) 


CDF 


m 




Ivqq 


1.1 






20.2 ± 4.5(stat. ± syst.) 


DO 


[73] 


PP^VV + X {V ^W, Z) 


ivqq, vvqq 


3.5 


16.8 ±0.5 


[64] 


18.0 ± 2.8(stat.) ± 2.4(syst.) ± l.l(lumi.) 


CDF 


[8^ 



TABLE XX A tabular summary of EW gauge boson cross section measurements and theory expectations during Run II and 
the Fermilab Tevatron. 



Mode 


TGC(s) 


jLdt A 
(fb-i) (TeV) 




95% CL Limits 

Ak-j 


1 (model constraints in parentheses) 
Az Agf Akz 




Expt. 
Ref. 


Ivy 


WW-y 


~ 0.7 


2.0 


-0.12,0. 


13 


-0.51,0.51 










DO [48] 


ivlv 


ww{n,z) 


3.6 


1.5 


-0.16,0. 


16 


-0.63,0.72 




-0.24,0.34 


Agi - Ak-, tan^ t 


■)w 


CDF 








2.0 


-0.14,0. 


,15 


-0.57,0.65 


A'y 


-0.22,0.30 


Agi - Ak-, tan^ t 


■)w 


CDF 






~ 1.1 


2.0 


-0.14,0. 


,18 


-0.54,0.83 


A^ 


-0.14,0.30 


Agf - AK-y tan^ t 


■)w 


DO [63] 








2.0 


-0.14,0. 


,18 


-0.12,0.35 


A^ 





AKj 




DO [63] 


mv 


wwz 


~ 1.0 


1.5 








-0.18,0.22 


-0.15,0.35 


-0.14,0.31 






















{Akz = Afff ) 




DO [67] 








2.0 








-0.17,0.21 


-0.14,0.34 


-0.12,0.29 
























{Akz = Apf ) 




DO [67] 


Ivqq, ivy 


WW{f,Z) 


0.4 


1.5 


-0.18,0. 


,17 


-0.46,0.39 


A-y 





AKj 




CDF [75] 






1.1 


2.0 


-0.10,0. 


,11 


-0.44,0.55 


A-y 


-0.12,0.20 


Agi - Ak-, tan^ t 


■)w 


DO [76] 










-0.11,0. 


,11 


-0.16,0.23 


A^ 





AKj 




DO [76] 


Ivy, lviv,mv, ivqq WWi^l, Z) 


~1.0 


2.0 


-0.08,0. 


,08 


-0.29,0.38 


A^ 


-0.07,0.16 


Agi - Ak-, tan^ i 


■)w 


DO [82] 










-0.08,0. 


,08 


-0.11,0.18 


A^ 





AKj 




DO [82] 



TABLE XXI A summary of anomalous charged TGC limits for different diboson processes and decay channels measured during 
Run II and the Fermilab Tevatron. 



CDF and DO. As evidenced from Fig. |31 constraints on 
the Higgs mass are limited more by the uncertainty in the 
W boson mass than by the uncertainty in the top-quark 
mass, which reflects the tremendous progress than has 
been made on the top-quark mass measurement in Run 
II. In fact, the top-quark mass is now the most accurately 
known quark mass (as a percentage of its value) , surpass- 
ing even its SU (2)l partner the h quark, whose mass is 
known to about 2.6% {mb{ihb) = 4.24 ± 0.11 GeV/c) 

The top-quark mass measured at the Tevatron corre- 
sponds closely to the pole mass. The pole mass is the 
mass the quark would have in the absence of confinement. 
Although the top quark decays on a time scale less than 



the time scale associated with confinement (Ag^^), the 
top-quark pole mass is nevertheless affected by confine- 
ment 84] . The ambiguity in the top-quark pole mass is 
of order Aqcd ~ 200 MeV, which is considerably less 
than the present uncertainty. A more precise definition 
of mass will be needed if one attempts to measure the 
mass to an accuracy much less than 1 GeV/c^. 

Both CDF and DO have made measurements of the top 
quark mass in a variety of final states and using differ- 
ent methods. The general procedure is to select events 
consistent with it production in either the dilepton or 
lepton plus jets final state. In the SM, top quarks decay 
via t -> Wb with a branching fraction of essentially 100%. 
The dilepton final state is that in which both W bosons 
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Mode 


TGC(s) 


jLdt 
(fb-i) 


A 
(TeV) 




95% CL Limits 




Expt. 
Ref. 


mi 


ZZZ, ZZ-^ 


1.0 


1.2 


-0.28 < /2o < 0.28 


-0.26 < /2o < 0.26 -0.31 < f^o < 0.29 


-0.30 < /5I, < 0.28 


DO f69j 






1.0, 3.6 


1.5 


\hJo\ < 0.033 


IftJol < 0.0017 |/ifo|< 0.033 


Iftfol < 0.0017 


DO [58] 



TABLE XXII A summary of anomalous neutral TGC limits for different diboson processes and decay channels measured during 
Run II and the Fermilab Tevatron. 
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FIG. 38 The 68% CL and 95% CL two-parameter limits on 
the ^WW/ZWW coupling parameters Ak^,, A, and Agf , in 
the LEP parametrization scenario and Anp ~ 2 TeV for the 
DO analysis [t^. The dots indicate the most probable val- 
ues of anomalous couplings from the two-parameter combined 
(electron-|-muon) fit and the star markers denote the SM pre- 
diction. 
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FIG. 39 The 68% CL and 95% CL two-parameter limits on 
the 'yWW/ZWW coupling parameters Ak and A, in the equal 
couplings scenario and Anp = 2 TeV for the DO analysis 
[7^ . The dot indicates the most probable values of anomalous 
couplings from the two-parameter combined (electron-fmuon) 
fit and the star marker denotes the SM prediction. 
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FIG. 40 Data compared with the sum of the predicted EW 
and multijet backgrounds for the Ai/i-lS* variable in the CDF 

analysis 81]. The band represents the total systematic uncer- 
tainty on the background. The measured signal is included 
here in the EW contribution. 



from the ti pair decay to leptons, and the lepton-f-jets 
final state is that in which one W decays to leptons and 
the other to quarks. Measurements of the mass have also 
been made in the all jets final state in which both Ws de- 
cay to quarks, although this offers less precision because 
of significantly higher backgrounds. 
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FIG. 41 Top: Comparison between data and fitted back- 
ground only. The measured signal is shown unstacked. The 
band represents the systematic uncertainty due to the shape 
of EW background as described in [sj]. Bottom: Compar- 
ison of the diboson signal (solid line) with the background- 
subtracted data (points). The dashed lines represent the ilcr 
statistical variations on the extracted signal. The gray band 
represents the systematic uncertainty due the EW shape. 
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FIG. 42 A graphical summary of EW gauge boson cross sec- 
tion measurements and theory expectations during Run II 
and the Fermilab Tevatron. 



An initial selection of candidate dilepton events is 
made, typically requiring two high px charged leptons, 
significant If^j^ attributed to neutrinos, and at least two 
jets. Typical thresholds for the lepton and jet transverse 
momenta range between 15 GeV/c and 20 GeV/c. The 
If^rp is attributed to the two neutrinos from ti decay, and 
the thresholds are higher, ranging between 35 GeV and 



50 GeV depending on topology. The initial selection of 
lepton-|-jets candidate events requires a single charged 
lepton, $rp and at least three (and often four) jets. The 
lepton and jet thresholds are similar to the dilepton final 
state, but the Iflj, threshold is typically relaxed giving re- 
quirements ranging from Ipj. > 15 GeV to $rp > 20 GeV. 
Because the lepton-|-jets events have significant back- 
ground from VF-|-jet events, an additional requirement 
that one or more of the jets is consistent with B hadron 
production (^-tagged) is often made. The typical tag- 
ging efficiency is 50%/jet with light quark misidentifica- 
tion rates of roughly 1%. The probability to misidentify 
a charm quark induced jet as being consistent with B 
hadron production is roughly 15%. Actual values vary 
somewhat depending on the 6-tagging algorithm used. 
The values also depend on the signal-to-background re- 
quirements of the specific analysis which allow more or 
less restrictive tagging requirements. These loose initial 
selections have minimal top quark mass bias and form 
the basis for further analysis. 

Typical signal and background estimates are given in 
Tab. lXXIIII for various channels. In the selections, trade- 
offs are made between signal acceptance and purity. In 
particular, the requirement that an event have one or 
more jets consistent with being initiated by 6-quarks con- 
siderably improves the signal purity. 

After the initial selections, various methods are used to 
determine the top mass. Among these are: (1) template 
methods in which data are fit to the predicted signal plus 
background distribution for a variable that is sensitive to 
the top mass, (2) ME methods similar to those discussed 
in Sec. lIII.E.4l in which tt event kinematics determined by 
matrix element calculations are convoluted with detector 
resolution functions to predict measured kinematic distri- 
butions as a function of the top mass used in the matrix 
element calculations, and (3) other weighting methods in 
which missing information is supplied by trying possible 
reconstructed jet and true ti decay parton assignments 
and computing a probability for the particular assign- 
ment as a function of top mass. The ME and various 
weighting methods compute event- by- event likelihoods 
as a function of top mass, and then compute the total 
likelihood for a data set as the products of the event like- 
hhoods. Combinations of these methods are also used. 

For lepton-|-jets events, the momenta of the four jets 
and the lepton are fully measured, and the neutrino mo- 
mentum in the plane transverse to the beam is inferred 
from the event . This gives seventeen measured values. 
Two constraints are available. The first is that the invari- 
ant mass of one pair of reconstructed jets should be con- 
sistent with that of the W boson, and the second is that 
the reconstructed mass of the two top quarks should be 
consistent with each other. This information allows a one 
degree-of-freedom fit to be performed in which the top 
mass is a free parameter. The resulting event top mass 
values are not simply averaged to determine the mass, 
but provide an input variable to the various methods de- 
scribed above. The fit can also be used to determine 
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JjCdt 


Expected 


Estimated 


Observed 


fe-tag 


Expt. 


Channel 


(fb-i) 


Signal Yield 


Total Yield 


Yield 


Required? 


Rcf. 


U 


2 


43.8 ±4.4 


215.8 ± 21.9 


246 


= 


CDF[85J 






78.0 ±6.2 


97.5 ±7.2 


98 


> 1 




££, efj. 


1 


36.7 ±2.4 


44.5 ± 2.7 


39 


no 


D0[86] 


ee 




11.5 ± 1.4 


14.8 ± 1.5 


17 










8.3 ±0.5 


13.7 ±0.7 


13 






e+track 




9.4 ± 0.1 


10.3 ± 0.2 


8 


> 1 




/i+track 




4.6 ±0.1 


5.5 ±0.1 


6 






£ + jets 


1.9 


183 ± 25 


247 ± 29 


284 


= 1 


CDF [87] 






69± 11 


75± 11 


87 


> 2 




£ + jets 


1 


162 ± 11 


57.2 ±4.2(est) 


200 


> 1 


D0[88J 



TABLE XXIII Signal and background yields for preselected tt samples from a subset of the published results. The purity 
varies considerably depending on the final state and whether or not fe-tagged jets are required. 



(or weight) the most likely parton-jet assignment in an 
event. Complications to this fit arise from the presence 
of intrinsic transverse momentum of the tt system and 
from final state gluon radiation giving rise to additional 
jets. The can be incorporated by adding an addi- 
tional variable for the unknown transverse momentum X 
to the fit and constraining it using a priori estimates of 
the spectrum. The impact of final state radiation can be 
controlled by requiring exactly four reconstructed (good) 
jets in the event. 

For dilepton events, the number of measurements mi- 
nus the number of constraints does not give enough infor- 
mation to fully determine the final state kinematics, and 
one assumption must be made. A variety of methods 
are used including ME weighting and neutrino weight- 
ing. The ME weighting method [89l - [9]| is related to the 
general ME method described below (and also used for 
^+jets analyses) . The neutrino weighting method [o^; Isl] 
is unique to dilepton events. In this method, an event 
weight is defined as a function of hypothesized top quark 
mass Mt using a comparison of the measured l^j, to ]^rp. 
values predicted for a set of possible neutrino pseudo- 
rapidity values 77j- , 77^^^ . Large weights correspond to 

situations in which the measured and predicted frp are 
similar and thus give a probability for different rj values 
for the neutrinos in each event. 

The use of the ME approach has become widespread, 
especially for £+jets final states. It is used either as the 
final analysis variable used to determine the measured 
mass or to provide information used in joint likelihoods 
for determining the top mass. In this approach, the prob- 
ability that the the jet momenta, lepton momentum and 
missing transverse energy observed in a given event, as- 
suming it arises from it production and decay, is com- 
puted bvjssj 

^« = ^ / E '^^(y; ^^t) /(gi) /(gs) W{y, x) 

in which da is the differential cross section for production 
of the final state partons in (momentum) configuration 



y for a given Mt, /(gi) and /(g2) are the parton density 
functions for the proton and antiproton, W{y, x) is the 
transfer function for computing the probability that a 
partonic final state y gives the observed reconstructed fi- 
nal state X and is a normalization factor equal to the 
expected observed cross section for a given Mt . For each 
event this probability, and the corresponding probabil- 
ity computed assuming the event is a background event, 
Pbkgix), are combined to give an event probability as a 
function of top mass 

P{x; Mt) - A{x)ifPt-t{x; Mt) ± (1 - f)Pbkg{x)) 

in which A(x) is a normalization factor incorporating ac- 
ceptance and efficiency effects and / is the fraction of ti 
events in the total sample. The overall normalization of 
P is forced to unity so that it can be used as a probabil- 
ity density. The event probability is then used either as 
input to additional likelihoods or the joint likelihood for 
the entire sample is formed as the product of the event 
probabilities and minimized as a function of Mt- 

The methods used in early Run II top mass measure- 
ments had systematic uncertainty from the energy cali- 
bration of reconstructed jets which, if unmodified, would 
have quickly surpassed the statistical uncertainty. A 
method was developed to incorporate an in situ calibra- 
tion scale factor as a second fit parameter when deter- 
mining the top mass. As an example, the ME method 
event probability definition was modified to the form now 
in general use. This form is 

P{x;Mt,kjes) = A{x) 

X [fP,j{x; Mt, k,es) ± (1 - f)Pbkaix; fc,e,)] 

in which kjes is a scale factor applied to all jet energies. 
In this form the transfer functions also become depen- 
dent on the calibration, W{y,x,kjes)- The negative log 
likelihood minimized to determine the top mass thus be- 
comes a function of two variables, mtop and kjets, both of 
which are varied during the minimization. As an exam- 
ple, Fig. shows the 2D contour in the Aj^s = kj^s — 1 
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FIG. 43 Contours of constant -log(L) in the A^es = kjes — 1 
versus Mt plane[93]. A perfect a priori jet energy calibra- 
tion used as input to this analysis would result in kj^s ~ 1 
corresponding to Ajes = 0. 

versus mtop plane from the Mt measurement using the ££ 
and ^-t-jets final states simultaneously. 

Systematic uncertainties in top mass measurements 
arise from both experimental sources and theory. Al- 
though the contribution from calibration uncertainties 
has been significantly reduced as described in the pre- 
vious paragraph, the dominant contributions to the sys- 
tematic uncertainty from experimental sources are the 
absolute reconstructed jet energy calibration and the rel- 
ative calibrations of jets initiated by 6-quarks and those 
from light quarks (u, d, s). Of these, the uncertainty 
from the jet calibrations dominates the experimental un- 
certainties. Uncertainties arising from theory include 
production model as assessed by comparing event gener- 
ators and fragmentation. The production model uncer- 
tainty becomes the dominant theory uncertainty when 
various measurements are combined. 

Tab. IXXIVI and Fig. give a summary of the pub- 
lished top mass results and a comparison with the 2008 
Particle Data Group (PDG)[33 average. When multi- 
ple results from one experiment using the same method 
and final state were available, only the highest integrated 
luminosity measurement is reported. The result shown 
in Tab. IXXIVI with highest integrated luminosity used a 
J C dt — 2.9 fb^^ data sample. The Tevatron experi- 
ments have now reported preliminary results using up to 
4 fb~^, and these results have been included in a world 
average combination by the Tevatron Electroweak Work- 
ing Group, with the most recent such combination josj ] 
giving 

Mt = 173.1 ± 0.6 (stat.) ± l.l(syst.) GeV/c^ 

For this combination, which does include preliminary re- 
sults, the systematic uncertainty is significantly larger 
than the statistical uncertainty. The dominant compo- 
nents of the systematic uncertainty arise from uncer- 
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FIG. 44 Published top mass measurements. For a given ex- 
periment and final state, the results from different methods 
have a significant correlation. Also shown are the most recent 
Tevatron combination 95] and the PDG world average. Both 
of these include a subset of the measurements listed in this 
table. 



tainty on the in situ jet energy calibration, uncertainty 
related to the tt production and decay model and uncer- 
tainty arising for color recombination. These contribute 
0.48 GeV/c^ 0.49 GeV/c^ and 0.41 GeV/c^ to the un- 
certainty, respectively. 

An alternative method of measuring the top quark 
mass makes use of the dependence of the production cross 
section on the mass. The top quark MS mass may be 
extracted directly from such a method [96j . 

By the end of the Tevatron running, the experiments 
expect to increase the data set by at least a factor of four 
compared to the published results and at least a factor of 
two for the data sets used in the combination. Thus, the 
statistical precision can be improved, but the systematic 
uncertainties will also need to be improved if the full 
statistical power of the final data set is to be realized. 



B. W Helicity from Top Quark Decays 

In the SM, top quarks decay with a V-A interaction 
through the process t W^b with B{t W + b) ^ I. 
The SM predicts a fraction /o = 0.697 ± 0.012 of top 
quark decays give a longitudinally polarized W and a 
fraction /+ = 3.6 x 10^** give a right-handed W . This 
leaves /_ — 1 — /q — /+ giving left-handed Ws. Gompared 
to the decay of lower mass quarks, the fraction of decays 
with longitudinally polarized W is significantly increased 
because of the large Yukawa coupling between the top 
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J C dt Analysis Main Systematic Uncertainty Expt. 



Channel 


(fb-i) 


Method 


Mt (GeV/c^) 


Source 


Value 
(GeV/c2) 


Ref 


U+{e + ]ets) 1.9 




171.9 ± 1.7(stat + JES) ± l.l(syst) 






CDF [94] 


U 


2.9 

2 

1.2 

1 

1 


(f) WT 
Matrix 

Template+Mt 
Matrix 

Matrix+z/WT 


165.5+^Kstat) ± 3.1(syst) 
171.2 ± 2.7(stat) ± 2.9(syst) 
170.71:1.9 (stat) ± 2.6(syst) ± 2.4(theo) 
164.5 ± 3.9(stat) ± 3.9(syst) 
174.7 ± 4.4(stat) ± 2.0(syst) 


jet calibration 
jet calibration 
jet calibration 
jet calibration 
jet calibration 


2.2 

2.5 

1.8? 

3.5 

1.2 


CDF [92 
CDF[85] 
CDF|98] 
CDF [99] 
DO [86] 


e + jets 


1.9 
1 

0.7 
0.4 


Matrix 
Matrix 

Decay Length 
Ideogram 


172.7 ± 1.8(stat) ± 1.2(syst) 
171.5 ± 1.8(stat + JES) ± l.l(syst) 

180.7li^-^(stat)±8.6(syst) 
173.7 ± 4.4(stat + JES) ± 2.1(syst) 


generator, jet calib 
b/u, d, s, frag. 

Bkg. shape 
b/u, d, s calib. 


0.6, 0.5 
0.83 
6.8 
1.7 


CDF [87] 
DO [88] 
CDF[100] 
D0[101] 




0.3 


Template+DLL 


173.5t^;^(stat + JES) ± 1.3(syst) 


signal model 


1.1 


CDF[102] 


AH jets 


1 

1.0 


Lineshape 
Matrix 


174.0 ± 2.2(stat) ± 4.8(syst) 
171.1 ± 3.7(stat + JES) ± 2.1(syst) 


jot calibration 4.5 
parton shower, jet calib 0.6, 0.5 


CDF[103] 
CDF[104] 




0.3 


Spectra 


172.3"9*6(stat) ± 10.8(syst) 


jet calibration 


9.6 


CDF[105] 


Tevatron (inch prel." 




173.1 ± 0.6(stat) ± l.l(syst) 






CDF+D0[95] 


Particle Data Croup 2008 


171.2 ± 2.1 






PDG[37] 



TABLE XXIV Summary of published Run II top quark mass measurements. Additional preliminary results using up to 
J jC dt ^ 4 fb"^ have been reported by both CDF and DO. The Particle Data Group average includes only some of these 
results, and the Tevatron average is based on a subset of the results in this table and additional preliminary results. 



quark and Higgs boson. Changes to the Lorentz {V — A) 
structure of the decay will result in changed values for 
the polarizatio n fraction s. Indirect results from 6 — 37 
measurements |l06l : Il07l | constrain /+ to less than a few 
percent. 

The angular distribution of the electron-type (It, — 
— 1/2) W decay product in < — > Wb is given by 

^{c} - ^[2(1 - c')fo + (1 - cff- + (1 + cff+] 

in which c = cos 9* is the cosine of the decay angle in the 
W rest frame meas ured with res pect t o the top quark di- 
rection. Both CDF [Toi] and DO [lo^] have measured the 
helicity fractions in top quark decay. The CDF analysis 
uses J C dt ~ 1.9 fb~^, and the DO result is based on 
/£ dt = 1.0 fb-i. 

Both experiments select events in the i+jets topology, 
and DO also uses the ££ + X topology. The initial selec- 
tions are similar to those described earlier in Sec. IIV.AI 
CDF uses a second selection with somewhat different 
requirements to perform an independent analysis. The 
CDF results are from the combination of the two meth- 
ods. 

All results rely on reconstructing the cos 6'* distribu- 
tion and comparing that to a set of predicted distribu- 
tions each of which is generated with a different pair of /o 
and /-!_ values. Several ambiguities arise in reconstructing 
cos 9* . These include the assignment of the reconstructed 
jets to the quarks from W decay (for the hadronic side 
of £-|-jets events) and the impact from the unmeasured 77 
value for the neutrino arising in leptonic W decay. These 
issues are handled differently in the three (2 CDF, 1 DO) 
analyses. The CDF analyses use only the leptonic W de- 
cay. The DO analysis also uses the jets from the hadronic 



W, but only | cos 6**1 is reconstructed for these, not cos 6'*. 
Although this loses the ability to distinguish from /_ 
for the hadronic decay, /o is better constrained. 

The results are extracted using likelihood fits to test 
cos 61* distributions. For one of the CDF methods and 
the DO analysis, predicted reconstructed cos 9* distribu- 
tions including signal and background contributions are 
generated using simulated events. The distributions are 
generated for a range of /o and f+ values, and best values 
are extracted using likelihood comparisons to the data. 
For the second CDF method, the result is compared to a 
distribution generated by convoluting the true cos 9* dis- 
tribution divided into six bins with a migration function 
which gives the probability that an event generated in a 
given true bin ends up a given reconstructed bin. This 
function includes the effects of resolution and acceptance. 
Effects of helicity on the acceptance are considered in all 
cases. 

The experimental results are summarized in 
Tab. IXXVI The results are consistent with each 
other and with the SM prediction. They are are 
generally limited by the statistical precision of the data 
set and are expected to improve by a factor between 
x2 and xlO by the end of Run II, depending on the 
analysis and machine performance achieved. Even if 
there is no improvement in the systematic uncertainty, 
the final results will still be statistics limited, but the 
statistical and systematic uncertainties will be similar. 
The dominant systematic uncertainty in the DO result 
comes from the ti production and decay model which 
was tested by comparing res ults from events generated 
using the pythia generator |llC| aii d with the default 
events generated using ALPGEN |111| . The CDF result 
does not have a single dominant source of systematic 
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CDF 


DO 


/o 




0.66 ±0.16 ±0.05 


0.425 ±0.166 ±0.102 


/+ 




-0.03 ±0.06 ±0.03 


0.119 ± 0.090 ±0.053 


fo, f+ 


= 


0.62 ±0.10 ±0.05 


0.619 ±0.090 ±0.052 


/+ 1 /o = 


0.70 


-0.03 ±0.04 ±0.03 


-0.002 ± 0.047 ± 0.047 



TABLE XXV W boson helicity fractions in tt decay deter- 
mined from fits to data. The first two lines are the results 
from simultaneous fits for /o and /+. The third line is the re- 
sult from the fit for fo with /+ = 0, and the fourth is the from 
the fit for /+ with fo set to the SM value. In all cases, the sta- 
tistical uncertainty is shown first and the second uncertainty 
is the systematic uncertainty. 




FIG. 45 The s-channel (left) and t-channel (right) diagrams 
for single top production. 



uncertainty. Although the contributions vary somewhat, 
the uncertainties arising from jet calibration, the back- 
ground prediction and modeling of final state radiation 
arc roughly equivalent to each and provide most of the 
systematic uncertainty. 



C. Single Top Production and Vtb 

The top quark was first observed in reactions medi- 
ated by the strong interaction process pp — >■ ti [ll'j Ill3l | , 
and nearly all measurements of the top quark properties 
have been made using such events. However, t op quarks 
can also be pro duced si ngly via the s-channel 114l4ll6j | 
and i-channel (I16l4ll9] electroweak diagrams shown in 
Fig. 03 These production modes involve the Wtb cou- 
pling and therefore cross section measurements provide 
previously unmeasured information on the electroweak 
sector, particularly the CKM matrix element Vtb- In 
addition to providing information about Vtb, detecting 
these events demonstrates progress toward Higgs sensi- 
tivity for the WH associated production mode with the 
subsequent decay H bb. 

The NLO cross sections for these processes ar e Cs — 
1.12±0.05pb and o-t = 2.34±0.13pb 0; [llMli for the 
s-channel and i-channel reactions, respectively, assuming 
|Vtf,| = 1 and a top quark mass of Mt — 175 GeV/c^. 
While the total single top cross section is slightly more 
than 40% of the tt cross section, the reduced parton mul- 
tiplicity in the single top events implies a signal-to-noise 
which is roughly 25 times lower than in tt making these 



events much more difficult to identify. In the SM with 
three generations, existing measurements of \Vub\ and 
I Vet I tightly constrain |Vtf,| to almost exactly unity [s^j- 
If there are more than three generations, however, |Vtfc| 
is almost entirely unconstrained. 

Both DO and CDF have long histories of searching for 
single top production. The m ost recent results from both 
DO JMliO] and CDF [HH give clear evidence for single 
top production and include results for the cross section 
and the CKM element |Vtb|. The DO analysis uses a data 
sample corresponding to J C dt = 2.3 fb~^. The CDF 
analysis uses up to J C dt — 3.2 fb~^. 

Candidate single top events are selected in the lep- 
ton+jets channel in which the W decays via either W — 
ev oyW ^ [IV. The event topology is thus one high-py 
charged lepton, corresponding to the neutrino, and 
two or more jets. Events are required to have fired single 
electron, single muon, electron+jet or muon+jet triggers. 
Both experiments use events with 2 or 3 reconstructed 
jets, and DO also makes use of 4-jet events. Tab. I XXVII 
shows the yields for the DO and CDF initial selections in 
which events are required to have at least one jet having 
a 6-tag. 

The final results for both experiments are based on 
combining results from multivariate (MV) techniques 
which exploit correlations among variables and event 
weighting methods based on calculated matrix elements 
for signal and background events. The methods are ap- 
plied to events selected by the initial offline requirements. 
DO uses three methods: (1) a ME method, (2) a Bayesian 
Neural Network (BNN) method and (3) a boosted deci- 
sion tree (BDT) method, and the final result comes from 
the combination of the individual results including cor- 
relations. CDF uses five methods: (1) a ME method, 
(2) a joint likelihood (LF) method, (3) a neural network 
method, (4) a boosted decision tree (BDT) and (5) pro- 
jected likelihood functions (LFS), and the final result is 
the combination of these including correlations among 
them. 

Final results are determined using likelihood fits. Like- 
lihoods are created from the multivariate outputs for each 
lepton flavor, jet multiplicity and 6-tag multiplicity for a 
given multivariate method. The joint likelihood formed 
by the product of the individual likelihoods for a given 
MV method is then maximized as a function of the signal 
cross sections with systematic correlations between lep- 
ton species, jet multiplicity and 5- tag multiplicity taken 
into account. The BLUE method [25| is used by DO to 
combine the results from each MV method into a sin- 
gle result. T he C DF results are combined using a super 
discriminant [l3l| |. in this case an additional neural net- 
work trained on the outputs of the multivariate methods. 
Fig. l46l shows the output of the three multivariate classi- 
flers for DO, and Fig. HTI shows the same for three of the 
CDF classifiers. The cross sections for each MV method, 
the combined cross sections and the fit probabilities are 
given in Tab. lXXVIll and illustrated in Fig.|48l The CDF 
and DO results have also been combined by the Tevatron 
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DO 




CDF 




fCdt 




2.3 fb~^ 




3.2 fb"^ 


2.1 fb"^ 


Process 


— 2 jets 


— 3 jets 


> 4 jets 


^ + jets 




tb + tqb signal 


139 ± 18 


63± 10 


21 ± 5 






tb signal 








77.3 ± 11.2 


29.6 ±3.7 


tqb signal 








113.8 ±16.9 


34.5 ±6.1 


W + jets 


1,829 ± 161 


637 ± 61 


180 ± 18 


1551.0 ±472.3 


304.4 ± 115.5 


Z+jets, 


229 ± 38 


85± 17 


26 ±7 


52.1 ±8.0 


128.6 ± 53.7 


Diboson 


in ^+jets 


in Z-^-^Qts 


in Z-\-^Gts 


118.4 ± 12.2 


42.1 ± 6.7 


tt 


222 ± 35 


436 ± 66 


484 ± 71 


686.1 ±99.4 


184.5 ± 30.2 


Multijet 


196 ± 50 


73± 17 


30 ±6 


777.9 ± 103.7 


679.4 ±27.9 


Total Pred 


2615 ± 192 


1, 294 ± 107 


742 ± 80 


3376.5 ± 504.9 


1404 ± 172 


Yield 


2579 


1,216 


724 


3315 


1411 



TABLE XXVI The integrated luminosity, predicted signal, backgrounds and event yields in the W+jets(6-tag) samples used 
in the CDF and DO single top analyses. 
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FIG. 46 The combined discriminant for the DO analysis. 



Electroweak Working Group. [132j 

The SM cross section is proportional to |Vtbp, so the 
cross section results are used to measure \Vtb\- The only 
assumption is that \Vtb\ ^ |Vts|, |Vtd| such that there is 
not a significant contribution to the signal from initial- 
state s or d quarks being transformed into top quarks 
via a i-channel W boson. These results are shown in 
Tab. IXXVIIll The measured single top cross section is 
in good agreement with the SM prediction with |Vtb| — 1- 

The precision of these results is limited by the statis- 
tics of the data samples. Both experiments are expected 
to continue these analyses through the end of the Teva- 
tron running. In the most optimistic scenario this will 
result in analyzed data samples of up to 10 fb^^. As- 
suming the results in Table IXXVllIl scale as C dt, 

the uncertainty on |Vth| will be reduced approximately 
two fold. 




0.2 0.4 0.6 0.8 1 

Super Discriminant 

FIG. 47 The combined discriminant for the CDF analysis 

V. HIGGS BOSON 

A. Precision Electroweak Constraints 



As discussed in Sections III.AI and III.Bl the standard 
model, with the S't7(2)L (8) U{1)y symmetry sponta- 
neously broken via the vacuum-expectation value of a 
single Higgs-doublet field, is consistent with precision 
electroweak data. These data are so precise that they 
are sensitive to the mass of the Higgs boson at one loop, 
despite the fact that the Higgs mass enters only logarith- 
mically. Fig. shows the constraints on the Higgs-boson 
mass from both precision electroweak measurements (the 
blue band) and direct searches (the yellow areas). The 
precision data depend on the extrapolation of the fine- 
structure constant a from its measured value at low en- 
ergy up to high energy, which suffers from an uncertainty 
associated with the contribution of low-energy QCD to 
the extrapolation. The solid black line indicates the cen- 
tral value, and the blue band takes into account all un- 
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rv „i — rr „ ~\~ rr^ 

^ st — s \ ^ t 




DO 


CDF 


ME 


4 40+^711 

^•^'-'—0.79 


^•^—0.8 


LF 






(B)NN 


4.70lJ-^^ 


1.8 ±0.6 


BDT 


q 7/1+0.95 


^•^-0.8 


LPS 




1 c+0.9 
J^-^-0.8 








Combined 


3.94 ± 0.88 


9 q+U.(i 
^•"^-0.5 


Predicted fl27] 


Tst = 3.46 


± 0.18 for Mt = 170 GeV/c^ 



[80] (j^t = 3.14 ±0.31 



TABLE XXVII Single top cross section measurements from 
the Tevatron (in pb). The DO result used Mt = 170 GeV/c^ 
in simulations and calculation. The CDP result used Mt — 
175 GeV/c. The predicted cross section changes by 0.1 pb 
for a 5 GeV/c^ change in Mt in this region. 




ME, DO 
BNN, DO 
BDT, DO 
Combined, DO [121] 

ME, CDF 
LF, CDF 
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BDT, CDF 
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I • iMET+jetS, CDF 

I Combined, CDF [122] 
Combined, CDF+DO [123] 
Theory [70; 112-119] 

I \ \ , \ \ \ l^*^^"?'^"'! I 

D123456789 10 

PIG. 48 Ideogram of the single top cross section measure- 
ments from DO and CDP. 



certainties. An alternative central value, associated with 
a different treatment of low-energy QCD, is indicated by 
the red dashed line. The effect on the central value by 
including the NuTeV data, which has some tension with 
the other precision electroweak data, is shown by the dot- 
ted magenta line. 



DO 


CDF 


IVt^l > 0.78 


\Vtb\> 0.71 


\Vtb\ = 1.07 ±0.12 


=0.9 ±0.11 



TABLE XXVIII \Vtb\ values extracted from the single top 
cross section measurements. The first row is the 95% CL 
result when no constraint is placed on on | Vti,|, and the second 
row is the result with the SM constraint < \Vtb\ < 1- 




itIh [GeV] 



PIG. 49 The A^^ as a function of assumed Higgs mass result- 
ing from fits to data assuming the SM. The exclusion regions 
from direct searches are also shown. 



It is striking that the precision data strongly prefer 
a Higgs boson mass in the 100 GeV/c^ region, while a 
purely theoretical analysis would acco mmod ate a Higgs 
boson as heavy as about 700 GeV/c^ [l33j |. Taken at 
face value, the precision data indicate the Higgs boson 
is not much heavier than the current lower bound of 
TJiH > 114 GeV/c^. It is also remarkable that the Teva- 
tron experiments have succeeded in excluding the region 
163 GeV/c^ <mH < 166 GeV/c^ 



B. Direct Searches 

The experimental strategy and sensitivity to discover 
the standard model Higgs boson at the Tevatron through 
direct searches depends strongly on the value of mn- 
For direct searches, the sensitive range is limited to 
mH ^ 250 GeV/c^ where there exists sufficient center-of- 
mass collision energy at the parton-level for an apprecia- 
ble rate of Higgs boson production via the gluon fusion 
(ggH) process which proceeds via a virtual quark loop 
(dominated by the top quark) (l34| . 

For niH 135 GeV/c^ (the "low mass Higgs" region), 
the dominant Higgs decay is to 6-quark pairs. The next 
largest contribution to the Higgs decay in the low mass 
Higgs region is —> rr which has a branching ratio of 
~ 7%. In the case of if — >■ bb, the ggH production is not 
useful to the direct Higgs search even using the bb invari- 
ant mass because of overwhelming continuum bb produc- 
tion. Instead, the associated vector boson (VH) process 
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where the Higgs boson is produced in association with 
a, W 01 Z boson is used. The additional leptons which 
can result from decay of the vector boson via W — )■ Iv, 
Z ^ U 01 Z ^ w with £ ^ e, fi [Hi; [l3i|, whether 
detected directly or indirectly through ^j,, can be used 
suppress backgrounds at the expense of a lower cross sec- 
tion for VH as compared to ggH. 

The cross section and cross section times branching 
fraction for s everal Higgs masses in the low mass Higgs 
region jl37| are shown in Tab. IXXIXI In addition to 
the inodes listed in Tab. IXXIXI the WH WWW 
£^£^ £~ channel provides some sensitivity in the transi- 
tion region around ttih ~ 135 GeV/c^ and vector boson 
fusion (VBF) qq — )■ qqH and ggH modes provide sensi- 
tivity in rf final states. 

For the low mass associated production channels VH 
in which the W or Z boson decay involves one or two 
charged lepton(s), an initial sample consistent with pro- 
duction of a vector boson decaying to charged and/or 
neutral leptons and having at least two jets is selected. 
This sample has very low signal-to-noise but is already 
dominated by vector boson events and therefore provides 
a control sample used to validate detector modeling and 
rate predictions of the dominant backgrounds in the final 
analysis. 

The requirement that one or more of the jets in the 
event is consistent with production of a B hadron is 
added to the selection, dramatically reducing the con- 
tribution from the initially dominant Vjj background 
and increasing the signal purity. Non-linear multivari- 
ate techniques for each final state are used to combine 
kinematic properties into a (typically) single variable in 
which a potential Higgs signal and backgrounds populate 
different regions. The distributions of these variables are 
then used as input to binned likelihood fits in which sig- 
nal and background fractions are allowed to fioat within 
constraints of the total yield, predicted background con- 
tributions and all uncertainties. The absence (presence) 
of a significant signal resulting from the fits determines 
the mass limit (indicates discovery). 

The low mass Higgs searches using the ZH — ?> i/i/bb 
channel and final states involving one or two r lepton de- 
cays follow different strategies. The ZH — )■ i/iybb channel 
has significant background from multijet QCD processes 
without real vector boson decay^ in which mismeasure- 
ment of jets results in significant In this channel, an 
additional step (described later in this section) is added 
to improve an understanding of this background. The 
search using final states that involve either one or two 
T leptons does not make use of 6-jet identification since 
the T leptons can arise either from Higgs boson decay or 
from W ov Z decay. 

The latest low mass Higgs searches from CDF and DO 



^ Such backgrounds are also present in the other channels, but as 
a much lower fraction of the total background. 



are described in Section IV.B.5I The most recent com- 
bination of low mass results, which includes some pre- 
liminary results not presented here, is described in Sec- 
tions ESI- EUl 

For friH ^ 135 GeV/c^ (the "high mass Higgs" re- 
gion), the dominant Higgs decay is to W boson pairs 
with the next largest contribution being 
Aside from a region around ttih = 150 GeV/c^ where 
BR{H ZZ'^*^) peaks to - 10%, the H WW'^*^ de- 
cay where one of the final state W bosons is off-shell for 
Mh less than twice the W mass completely dominates 
over the sensitive region of the Tevatron. The Higgs pro- 
duction cross section and branching fractions for several 
Higgs masses in the high mass region are shown in Ta- 
ble IXXXI The standard model Higgs boson branching 
fraction to WW* varies from 7.5% at 115 GeV/c^ to 
73.5% at 200 GeV/c^ with a maximum of 96.5% at sa 
170 GcV/c^ [133. For M// - 170 GeV/c^ the standard 
model Higgs boson decays almost exclusively to two on- 
shell W bosons, making this a region where the Tevatron 
has the best chance at a discovery or an exclusion. 

At the Tevatron, the most sensitive Higgs search chan- 
nel over the range 135 GeV/c^ ^ mn ^ 200 GeV/c^ is 
gg ^ H ^ WW'^*^ — >■ iviiy, where the two charged lep- 
tons in the final state are of opposite charge [138; 13.9j. 
The iiyii^ final state represents 6.0% of all WW* de- 
cays, where £ is either an electron or a muon, including 
those from r leptons produced in the W decays. Over 
the last few years, the experimental search for a high 
mass Higgs in decay to WW at the Tevatron has evolved 
substantially to include powerful multivariate techniques 
and additional H WW^*^ processes such as VBF 
H WW'^*^ and VH{^ WW'^*^). The latest high mass 
Higgs searches in the H — > Tyiy^*) channel from CDF 
and DO are described in Section IV.B.5I A combination 
of these these results are described in Section IV. CI 

The following sections describe the individual chan- 
nels used in the Higgs search starting with the low mass 
channels, WH £M, ZH Ubb, ZH wbb 
and T final states, and finishing with the high mass 
H WW'^*^ — >■ ivlv search. The Higgs boson section 
concludes with a presentation of the most recent set of 
CDF and DO combined limits. 



1. WH luhh Final State 

This final state has the largest cross-section times 
branching ratio of the entries shown in Tab. IXXIXI and 
it also gives the most sensitivity to Higgs boson pro- 
duction for Mh ^ 135 GeV/c^ Both CDF [m^ and 
DO fliU have carried out searches in this final state. 
The currently published results for DO use a sample of 
J£ dt ^ 1.0 fb-i. The CDF pubhshed results use a 
sample of J £ dt = 2.1 fb""'^, and preliminary results 
have been reported by both collaborations which use up 
to approximately J £ dt — 4 fb~^. The search strategies 
used by the two collaborations in the WH channel are 
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Tif'H l^vjG V / C ) 


G ypp — r VV 11 J (J X JD\^VV - — r j 




(T X JDyZj — T tZj (7 X iI>yZj — r UV j 


100 


286 63.6 


167 


11.2 33.4 


110 


209 46.4 


124 


8.35 24.8 


120 


153 34.0 


92.7 


6.24 18.5 


130 


115 25.5 


70.9 


4.77 14.2 


TABLE XXIX The production cross section and cross section times branching fraction 


for the low mass {niH < 135 GeV/c^) 


Higgs boson search. 


Cross sections are in femtobarns. For the final states denoted and 11, the branching fraction used is the 


sum of £ — e jj,. 








ra„ (GeV/c^) 


<^{99 ^ H) 


cr(VBF) 


B{H^WW) 


150 


548 


45.7 


0.682 


160 


439 


38.6 


0.901 


170 


349 


33.6 


0.965 


180 


283 


28.6 


0.934 



TABLE XXX The production cross section and cross section times branching fraction for the high mass (m/j > 135 GeV/c^) 
Higgs boson search. Cross sections are in femtobarns. 



generally similar. 

The initial event selection requires one high-pr elec- 
tron or muon, large ]^rp and two or more jets. This selec- 
tion mirrors the £, v and hh pair present in the final state. 
CDF requires the events to be selected by a trigger based 
on the presence of a high-energy electron or muon.'^ The 
DO analysis requires l — e events to be selected by single- 
electron or electron plus jet triggers. For the muon final 
state DO uses a two phase trigger selection. In the first 
phase, the events are required to be selected based on at 
least one of a set of single lepton or single-lepton plus 
jet triggers, and all yields are predicted and compared 
with data. After establishing agreement in this pass, the 
analysis is repeated allowing events to be selected by any 
trigger. The yield increases by an amount predicted by 
the inefficiency of the triggers for the first pass and all 
kinematic distributions remain in agreement after a sim- 
ple yield scaling to 100% trigger efficiency. 

At this stage, the dominant event source is VF-l-jets 
with the jets arising from light quark (u, c?, s, c) pro- 
duction. The signal-to-background is then improved by 
requiring one or two of the jets in the event to be iden- 
tified as consistent with h qua rk production. The DO 
h identification algorithm jl42| uses a Neural Network 
(NN). The NN is trained and verified on a combina- 
tion of data and simulation, with the critical efficiencies 
and misidentification rates determined from data control 
samples. The CDF analysis used three different h iden- 
tification algorithms. Jets are 5-tagged by one or more 
algorithms bas ed on the presence of a secondary vertex 
(SECVTx) [l43l | , a neu ral network [l44| , or signed impact 
parameters fJP) [l45l |. 



3 The event selection in recent preliminary results from CDF in- 
clude additional events in this channel which are selected by 
triggers. 



The signal purity is greatest when two h jets are re- 
quired, but this introduces significant efficiency loss com- 
pared to the case in which only one h jet is required. To 
gain the most sensitivity, both categories of events, single 
6-tagged and double 6-tagged, are retained, but they are 
analyzed separately. Optimization studies indicate that 
for the single-tag channel, rather restrictive tagging re- 
quirement is needed to control the background contribu- 
tion from events with light-flavor (If) jets misidentified 
as h jets while for the double-tag channel, a less strin- 
gent requirement suffices thereby giving increased per- 
jet tagging efficiency. The DO single-tag analysis uses 
a &-tag operating point which gives a typical efficiency 
of 48% with misidentification rate of 0.5% and an op- 
erating point for the double-tag analysis which gives an 
efficiency of 59% and a misidentification rate of 1.7%. 
The CDF event selection divides events into three ex- 
clusive categories based on which algorithms identify a 
jet as 6-tagged. The first category contains events which 
have two jets tagged by the SECVTX algorithm. The sec- 
ond category contains events not selected into the first 
category which have one jet identified by the SECVTX al- 
gorithm and one by the jet probability algorithm. The 
third category contains events which are not selected into 
either of the other two categories and have one jet iden- 
tified by both the SECVTX and NN algorithms. 

The sensitivity of the analyses are further improved 
by using kinematic properties of the events to distin- 
guish signal and background events. The most impor- 
tant variable is the mass of the dijet system correspond- 
ing to the Higgs decay. For single-tagged events, the 
mass is computed using the tagged jet and the highest 
jet remaining, and for the double-tagged events the 
mass is computed using the two tagged jets. In addition 
to the mass, other variables can also distinguish signal 
and background on a statistical basis. To make best use 
of these, both experiments use a neural network to en- 
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hance the separation of signal and background in the 
final step of the analysis. The CDF neural network has 
six input variables and the DO network have seven input 
variables.^ The variables used by each experiment are 
listed in Tab. IXXXII for comparison. Because of limited 
statistics in the M^+three-jet sample, DO uses the NN 
only for W+2 jet case and uses the dijet mass for the 
W+3 jet case. Backgrounds to the WH search include 
W+jets production (including jets arising from heavy fla- 
vor (hf) production), ti, single top, diboson production 
{WW, WZ and ZZ) and a small contribution from mul- 
tijet events in which either the lepton is actually a jet 
misidentified as a lepton or the lepton arises from heavy 
flavor decay and the jet energy is not reconstructed so the 
lepton appears to be isolated. Both experiments estimate 
the ti, single top and diboson yields using the product 
of the theoretical (N)NLO cross-section for each process, 
the luminosity and acceptance times efficiency for each 
process. Corrections are applied based on comparison of 
data and simulated control samples. The M^+jets back- 
ground cross sections are poorly known, and data-driven 
approaches are used to estimate these yields. The mul- 
tijet background is hard to model from simulation, so 
data-driven methods are also used for this background. 
The two experiments use different data-driven methods 
as outlined below. 

The total VF-|-jets background after 6-tagging in the 
CDF result is estimated separately for M^-l-lf events and 
H^+hf events. The VF-l-lf contribution is estimated by 
applying data-derived mistag probabilities to untagged 
M^-|-jets samples. Three different tagging algorithms are 
used in the analysis, and the mistag probability determi- 
nation differs for each of these. For the SEC VTX and JP 
tagging algorithms, the mistag probability |l43l : Il45| is 
derived using events with a negative decay length with 
a correction applied to account for the heavy-flavor con- 
tent of the control sample used to determine the mistag 
probability. For the JP algorithm, the mistag probabil- 
ity is parametrized as a function of 77, primary vertex z 
position, jet Et, scalar transverse energy, and vertex and 
track multiplicity. For the TV TV-based tagging algorithm, 
a light-flavor rejection factor derived from control data 
samples is used. The W^-|-hf contribution is determined 
by measuring the heavy flavor faction in Vl^-|-jets events 
and applying a 6-tagging efflciency to these events. The 
initial heavy flavor fraction is derived from ALPGEN and 
PYTHIA simulation and is corrected by a factor of 1.4±0.4 
derived from a jet control data sample. The 6-tagging 
efficiencies are determined from simulation and checked 
using control data samples. 

The W-|-jets background yields after 5-tagging in the 
DO result are fixed by normalizing the simulated events 



* Recent preliminary results has also used the matrix element 
method adapted from the top mass measurement and decision 
trees for the multivariate technique instead of neural networks. 



to the untagged VT-l-jets data after subtracting the other 
backgrounds from the data. The relative contributions 
of the VF-|-lf, W + cc, and W + bb in the untagged sample 
are fixed to the cross-section ratios predicted by mcfm. 
The yields in the 6-tagged sample are then computed by 
applying flavor-based (mis)identiflcation probabilities to 
the jets in simulated events. The probabilities (one for 
d- and s- quark initiated-jets, one for c-quark initi- 
ated jets and for for 6-qua rk in itiated jets) are derived 
using data control samples jl42| and are parametrized as 
functions of jet px and 77. 

The small background from multijet events is difficult 
to simulate accurately, so this component is also deter- 
mined from control data samples. For CDF, a control 
sample is selected using events which have non-isolated 
leptons and low and the yield in the signal sample is 
determined by extrapolating the yield from this sample 
into the signal region having isolated leptons and high 
$]rp. For DO, this background is determined by select- 
ing a multijet dominated control sample with kinematics 
similar to the WH events, and then applying a proba- 
bility that these events would be misidentified and ap- 
pear in the signal sample. The multijet control sample is 
selected by requiring lepton candidates which pass very 
loose isolation requirements, and the background yield 
is computed by applying the event-by-event probability 
that these loose-isolation events would pass the standard 
isolation requirement and thus appear in the signal sam- 
ple. 

Tab. lXXXIII shows the data yields and background and 
signal predictions for the CDF analysis, and Tab. lXXXIIII 
shows the yields for the DO analysis. The predicted Higgs 
yield includes not only WH events, but also a small con- 
tribution from ZH — > ££bb events in which one of the 
leptons {£) is not identified and thus generates fj^ and a 
final state consistent with the WH final state. 

Both CDF and DO derive their final result by com- 
paring the predicted and observed spectra using binned 
likelihoods. For both experiments, the single-tag and 
double-tag samples are separated and their likelihoods 
combined for the final results. DO further separates the 
samples into events with exactly two jets and events with 
three jets. When determining the final result, CDF uses 
the dijet mass as the input to their likelihood. DO uses 
the NN output for the two jet events and the dijet mass 
for three jet events. The likelihood methods are described 
in section rvTCI 

Systematic uncertainties are evaluated by both ex- 
periments for trigger efficiency, lepton identification effi- 
ciency, b~iet (mis) identification efficiency, the jet identifi- 
cation efficiency and jet energy calibration, multijet back- 
ground calculation method, luminosity, and theory cross 
sections used for background and signal event yield calcu- 
lations. CDF additionally reports uncertainties from par- 
ton density functions and initial and final state radiation 
modeling. DO reports an additional systematic from the 
M^-|-jets simulation derived by comparing shapes of dis- 
tributions of data and simulated events before 6-tagging. 
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CDF DO 



M,j,j+ - invariant mass of Ji and J2 


Mjj - invariant mass of Ji and J2 


and trie closest loose jet il iXHyJ, Jioose) < U.9 




ESt (loose jets) 


EriJi) 


Pt{Ji)+Pt(J2)+pT(£) - $j. 


Et{J2) 




Ai?(Ji, J2) 




A0(Ji, J2) 


Wrie) + Pt{$t) + Pt{Ji) + Pt{J2)\ 


|PtW+Pt(^t)I 



Ipt(Ji) + Pt{J2) 



TABLE XXXI The neural network inputs for the CDF and DO WH analyses. Here Ji denotes the i-th jet in a list ordered 
by jet Et in which Ji is the highest Et jet in the event. M]^™ is the mass of the lepton, and the jet (Ji or J2) which 
gives the lower mass value. Finally, Vmax is the three momentum of the neutrino in which \pz \ is the larger of the two values 
calculated when forcing the lepton and neutrino system to have a mass equal to the W boson mass. 







Central region 






Plug region 




Pretag Events 




32242 






5879 




6-tagging 


ST+ST 


ST+JP 


ST+NN 


ST+ST 


ST+JP 


ST+NN 


Mistag 


3.88±0.35 


11.73±0.92 


107.1±9.38 


1.00±0.18 


3.18+0.49 


28.47+3.30 


Wbb 


37.93±16.92 


31.15±14.03 


215.6±92.34 


7.40±3.96 


6.23+3.37 


43.09+12.33 


Wcc 


2.88±1.25 


7.87±3.43 


167.0±62.14 


0.96±0.49 


1.53+0.81 


33.37+9.55 


tt(6.7pb) 


19.05±2.92 


15.56±2.39 


60.68±9.30 


2.14±0.34 


1.79+0.31 


7.17+1.00 


Single top(s-ch) 


6.90±1.00 


5.14±0.75 


14.38±2.09 


0.69±0.10 


0.51+0.08 


1.53+0.20 


Single top(t-ch) 


1.60±0.23 


1.87±0.27 


29.57±4.33 


0.22±0.04 


0.24+0.04 


3.54+0.47 


WW 


0.17±0.02 


0.93±0.11 


15.45±1.91 


O.OliO.Ol 


0.12+0.04 


3.00+0.20 


WZ 


2.41±0.26 


1.84±0.20 


7.59±0.81 


0.58±0.06 


0.42+0.05 


1.62+0.09 


ZZ 


0.06±0.01 


0.08±0.01 


0.31±0.03 


O.OOiO.Ol 


0.01+0.01 


0.02+0.00 


Z -^TT 


0.25±0.04 


1.29±0.20 


7.27±1.12 


O.OOiO.Ol 


0.01+0.01 


0.24+0.03 


non-W QCD 


5.50±1.00 


9.55±1.73 


184.7±33.04 


1.16±0.44 


1.51+0.55 


18.34+5.54 


Total Background 


80.6±18.8 


87.0±18.0 


809.6±159.4 


14.2±4.0 


15.5+3.6 


140.4+16.9 


WH signal (120 GeV/c^) 


0.85±0.10 


0.60±0.07 


1.70±0.14 


0.09±0.01 


0.06+0.01 


0.20+0.01 


Observed Events 


83 


90 


805 


11 


13 


138 



TABLE XX XII The predicted and observed yields for the CDF 2 fb ^ WH search separated by detector region. This table 
is from Ref. [H^ 



W + 2 jet W + 2 jet W + 3 jet W + 3 jet 
16-tag 2 6-tag 1 6-tag 2 6-tag 



WH 


2.8 


+ 


0.3 


1.5 


+ 


0.2 


0.7 


+ 


0.1 


0.4 + 


0.1 


WZ 


34.5 


+ 


3.7 


5.3 


+ 


0.6 


9.1 


+ 


1.0 


1.7 + 


0.2 


Wbb 


268 


+ 


67 


54 


+ 


14 


87 


+ 


22 


22.7 + 


5.7 


ly+jets 


347 


+ 


87 


14.0 


+ 


4.4 


96 


+ 


24 


8.5 + 


2.7 


ti 


95 


+ 


17 


37.4 


+ 


7.0 


156 


+ 


29 


81 + 


15 


single t 


49.4 


+ 


9.0 


12.4 


+ 


2.3 


15.7 


+ 


2.9 


6.7 ± 


1.2 


m-jet 


104 


+ 


29 


8.9 


+ 


2.1 


54 


± 


15 


8.7 ± 


2.1 


Total 


896 


+ 


177 


132 


+ 


27 


418 


+ 


76 


129 + 


24 


Data 


885 






136 






385 






122 





TABLE XXXIII The predicted and observed yield s for the 
DO 1 fb'^ WH search. This table is from Ref. |l4ll | 



The CDF NN output distributions for the tagged event 
selections are shown in Fig.. [SO] The DO NN output 
distribution (two jet events) and dijet mass distribution 
(three jet events) are shown in Fig. [5T] The CDF and 
DO cross section limits for the WH final state only are 
shown if Figs. [52] and [53] respectively. The result from 
combining the results for all low mass final states is shown 



in Section IV^Cl 



2. ZH mb Final State 

The final state from ZH — ^ Mbb events has the lowest 
production rate among those shown in Tab. IXXIXI but 
has the highest purity. A similar strategy is taken for this 
final state as for the WH final state. An initial sample of 
Z+dijet events is selected, the purity is improved by re- 
quiring b jets, and the final result is determined by using 
either the invariant mass distribution or the output of a 
neural network as i nput s to a shape-based limit setting 
program. Both DO 146l| and CDF [li?] have published 
results in this final state. The DO publication uses a data 
sample corresponding to J C dt — 0.45 fb~^. The CDF 
published results corresponds to J £ dt = 1.0 fb^^. 

The initial selections for both experiments require 
events with two leptons whose invariant mass is con- 
sistent with a Z boson and two additional jets. The 
backgrounds in this channel arise from Z+jets, ti events, 
diboson production, and from jets misidentified as lep- 
tons or a lepton+jet system being misidentified as an 



42 



— • — Data 

I I W+HF 
I I Mistag 

tt (6.7pb),Singletop 

Diboson 
I I NonW 

Higgs(120GeV)x10 

^ Background Uncertainty 




0.7 0.8 0.9 
NN output 



Data 

W+HF 

Mistag 

tt (6.7pb),Single top 

Diboson 

NonW 

Higgs(120 GeV)x10 
Background Uncertainty 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

NN output 

FIG. 50 The NN output distributions for (a) single- 
(ST+NN) and (b) double-tagged (ST+ST, ST+JP) events 
from the CDF 2 fb"^ WH search. 



isolated lepton. They are determined using a combina- 
tion of simulated events and data control samples in es- 
sentially the same manner as for the WH channel. One 
difference with respect to the WH calculation is that 
sidebands of the ££ invariant mass distribution are used 
by DO to determine the misidentification and false isola- 
tion background. Both experiments identify b jets using 
secondary vertex algorithms^. The CDF result splits the 



^ For DO this differs from the WH case because the pubhshed ZH 
results predate the availabihty of the NN tagger. More recent 
DO preliminary results use the same NN tagger described in the 
WH analysis. 
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FIG. 51 (a) The NN output for double-tagged dijet events 
and (b) the dijet mass for double-tagged three jet events from 
the DO WH search. 



final tagged sample into single- and double-tagged sam- 
ples as in the WH analysis, but the DO result does not. 
The DO result is derived using the dijet mass distribu- 
tion only, while the CDF result is based on the binned 
likelihood of a NN output. The sources of systematic un- 
certainties for these channels are the same as for the WH 
result. Tab. IXXXIVI shows the yields for the single- and 
double-tag analyses from CDF, and Tab. IXXXVI shows 
the yields for the dijet and double-tagged samples for 
DO respectively. Figure [54l ((55)) shows the dijet invariant 
mass (neural network output) distribution from the DO 
(CDF) analysis. This distribution is used as the input to 
the limit setting program. 

Systematic uncertainties are included for trigger and 
lepton identification efficiencies, parton density func- 
tions, background cross sections, 6-tagging efficiencies, 
jet energy reconstruction and the methods used to es- 
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FIG. 52 The cross section limits from the CDF WH search. 




130 140 150 
Higgs Mass (GeV/c^) 



FIG. 53 The cross section hmits from the DO WH search. 



Source 


Single- Tagged 


Double- Tagged 


Z + bb 


35.1 ± 14.6 


6.3 ±2.5 


Z + CC 


21.8 ± 8.5 


1.0 ±0.4 


Z + qq 


32.3 ±5.5 


1.0 ±0.2 


tt 


5.2 ± 1.0 


2.8 ±0.6 


zz 


4.0 ±0.8 


1.3 ±0.3 


wz 


1.2 ±0.2 


0.04 ±0.01 


Non-Z 


1.9 ± 1.4 


0.2 ±0.2 


Total Expected 


101.5 ±32 


12.7 ±4.1 


Observed 


100 


11 


Signal Yield 


0.44 


0.23 


TABLE XXXIV Single- and double-tag 
CDF 1.0 fb"^ ZH ^ Ubb search. 


ged yields for the 



Source 


Dijet 


Double-tagged 


Z + bb 


17.4 


3.3 


z + jj 


851 


3.8 


tt 


12.3 


3.9 


wz + zz 


30.6 


0.74 


Non-Z 


44.1 


0.59 


Total Exp. 


956 


12.3 


Observed 


1008 


15 



TABLE XXXV Dijet and double-tagged sample yields for the 
DO 0.45 fb"^ ZH Ubb search. 
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FIG. 54 Dijet mass distribution for the DO 1.0 fb"^ ZH 
Ubb analysis. 
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FIG. 55 NN distributions for the (a) single-tagged and (b) 
double-tagged samples from the CDF 1.0 fb"^ ZH -5- Ubb 
analysis. The 95% CL upper bound cross section is 19 x the 
SM prediction. 
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FIG. 56 Expected and Observed limits for the CDF 1.0 fb~^ 
ZH — i> Uhh analysis. 



escapes detection. This is particularly true for the case 
I = [I because the muon leaves very little energy in the 
calorimeter and thus results in event similar to that 
from Z ^ vv decay. 

Unlike either of the previously discussed final states, 
the ZR vvhh final state has no charged leptons from 
vector boson decay. This implies a significantly increased 
background from SM multijet events in which arises 
from mismeasurement. This background is difficult to 
model from simulation, and analyses of this final state 
must develop techniq ues to measure it us ing data control 
samples. Both CDF and DO [Hil have published 

results in this final state. The two experiments have de- 
veloped different methods for controlling and estimating 
the multijet background. 
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FIG. 57 Expected and Observed limits for the DO 0.45 fb"^ 
ZH Mbb analysis. 



timate instrumental backgrounds. The CDF result also 
includes a systematic from the top mass uncertainty. The 
systematic uncertainties for the DO measurement range 
between 2% and 20% expressed as a fraction of the total 
background. The systematic uncertainties for the CDF 
measurement range between 1% and roughly 25% of the 
total background. The largest contribution for both ex- 
periments is from the background cross sections. 

The final limits, derived using the same methods as 
the WH results, are shown for CDF in Fig. and for 
DO in Fig. [ST] As expected these channels have consid- 
erably less sensitivity than the WH channels because of 
the significantly lower signal cross-section times branch- 
ing fraction. 



3. ZH 



ybb and Related Final States 



The ZH — >■ h'h'bb final state has a production rate in- 
termediate between the WH and ZH — > ££bb final states. 
This final state also has a significant contribution from 
the process WH — )■ (i/bb in which the charged lepton i 



a. CDF Search The CDF analysis uses a data sample 
corresponding to 1 fb~^ and begins with selection of 
events passing a l^j, trigger with level one > 25 
GeV, a level two requirement of two jet clusters having 
Et > 10 GeV and a level three requirement of J^rp > 35 
GeV. At least one of the level two jets must also satisfy 
rj < 1.1. The initial offline selection ("pre-tag") requires 
events to have MET > 50 GeV and exactly two jets with 
Et > 20 GeV. One of the jets must have Et > 35 GeV, 
and the other must have Et > 25 GeV. Additionally, 
one of the jets must satisfy \ri\ < 0.9, and the other 
jet must satisfy \r]\ < 2.4. The two jets must also have 
A(j) > 1.0 rad, and events with high pT, isolated leptons 
are vetoed. Finally, at least one of the jets is required to 
have a secondary vertex 6-tag 

All non-multijet backgrounds; ti, W+jets, Z-|-jets and 
diboson production are modeled using simulated events. 
The multijet background is studied by dividing the sam- 
ple into two control regions and a signal region. The 
regions are defined in Tab. IXXXVII The multijet back- 
ground in all regions is divided into two components: 
(1) events with only light-flavor (u, d and s) quarks in 
which one or more of the jets is misidentified with a sec- 
ondary vertex and (2) events with c and b quarks. The 
contribution from light flavor jets is determined using a 
control sample with no tags to which a misidentification 
factor is applied. The heavy fiavor contribution is deter- 
mined using simulated events with normalization factors 
for single-tag and double-tag topologies determined by 
forcing the data yields and the sum of all backgrounds 
to agree (before dividing the sample into the three re- 
gions). The scale factors are 1.30 ± 0.4 (1.47 ± 0.07) for 
the single(double) 6-tagged events. 

The final analysis selection is determined by optimizing 
S/ \/iB) where S is the total signal yield, including the 
contribution from WH — ^ ivbb in which the lepton is not 
identified, and B is the total background. The optimiza- 
tion is carried out only for the signal region. The final 
selection then requires A0(ji,j2) > 0.8, I^t/Ht > 0.454, 
> 60 GeV and ^t > 70 GeV. Here Ht is the scalar 
sum of the jet Et values, I^t is the magnitude of the 
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Region 


Dominant 
Source 


Selection 


CRl 


Muhijet 


Leptons Vetoed 
A${Et2,^t) < 0-4 


CR2 


EW, tt 


Lepton Required 
A^{Et2,$t) > 0.4 


Signal 




Leptons Vetoed 
A^{Eti,$t) > 0.4 
A^{Et2,$t) > 0-4 


TABLE XXXVI Definitions of the three regions (CRl, CR2, 

and Signal) in the CDF VH $j.bb search. A${Eti{2) , $t) 
is the angle between the Ifj. and the jet with the high- 
est(second highest) Et- 


Source 


Signal-tagged 


Double-tagged 


Muhijet 

tt 

Diboson 
W^+hf 
Z+hf 

Mistags 


93 ±23 

27.3 ±3.8 
7.0 ± 1.4 

33.4 ± 16.2 
18.3 ±8.1 

69 ±9 


3.74 ± 1.27 
4.88 ± 0.8 
0.79 ±0.19 
1.65 ±0.86 
1.67 ±0.77 
1.64 ±0.48 


Total 
Observed 


248 ± 43 
268 


14.4 ± 2.7 
16 



TABLE XXXVII Yields in the CDF VH $j.bb analysis. 



vector sum and Ji(j2) denotes the jet with the high- 
est(second highest) Et- 

The systematic uncertainties arise from a number of 
sources. For the CDF analysis the dominant uncertainty 
in the jet energy calibration which varies between 10% 
- 26% for multijet and F±jets (including heavy flavor) 
backgrounds, but is only 8% for signal events. The other 
dominant systematic arises from the calculated cross sec- 
tions used to normalize backgrounds. This ranges be- 
tween 11% - 40% for a given source depending on the 
samples. The total systematic is 17% for the single- 
tagged analysis and 19% for the double-tagged analysis. 

Fig. [SHI shows dijet invariant mass for the single- and 
double-tagged events in the signal region. The yields are 
given in Tab. IXXXVlH The limits on Higgs production 
are calculated using the same procedure as for the previ- 
ous two channels, and they are shown in Fig. 1591 



b. DO Search The published DO analysis uses a data 
sample corresponding to J £ dt = 5.2 fb~^. Events se- 
lected for this analysis must satisfy a trigger requiring 
two acoplanar jets and $!rp in which the kinematic thresh- 
olds varied as a function of the instantaneous luminosity. 
The primary requirements of the initial offline event pre- 
selection are at least two jets with Et > 20 GeV and 
\r]\ < 2.5, Acj){ji,j2) < 165°, and $t > 20 GeV. This 
is divided into four samples: (1) a signal sample defined 
by the additional requirements of $!t > 40 GeV, S* > 5, 
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FIG. 58 The dijet invariant mass distributions for the (a) 
single-tagged and (b) double-tagged CDF searches in the 
VH — )> l^j,bb final state. 



D < Tr/2 and a veto of events with high-p^, isolated lep- 
tons, (2) an electroweak test sample enhanced in W ^ fiiy 
events and defined similarly as for the signal sample, but 
requiring the presence of a high pT, isolated muon, (3) a 
multijet modeling sample defined similarly to the signal 
sample except D > Tr/2 and (4) a multijet enriched sam- 
ple used to confirm the background model predicted us- 
ing the multijet modeling sample and defined by 1/Jt > 30 
GeV and no requirement on S. Here S is the significance 

of the $Jt and D = A</i(.^-p, ^V) = 4"]/] ~ 4' J in which 

'frp is the magnitude of the vector sum of the transverse 
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FIG. 59 Expected and observed 95% CL upper bounds on 
the VH production cross section for the CDF l^j, + bb search. 



momentum of the tracks in the event and is the az- 
imuthal angle of either $!rp or "frp. 

The presence of b hadrons from the H bb signal 
decay is used to improve the signal to noise. Events 
were additionally required to have a least one of the two 
highest Et jets identified as being consistent with a b 
hadron jet as determined by a high purity (tight) 6-jet 
identification requirement. Events are then categorized 
as double-tagged if the remaining one of the two high- 
est Et jets satisfies a lower purity (loose) identification 
requirement. Events without an additional tag are de- 
noted single-tagged. For the tight requirement, the per 
jet identification efficiency is 50% with a misidentifica- 
tion probability of 0.5%. The loose requirement has an 
efficiency of 70% with a misidentification probability of 
6.5%. Tab . IXXXVIII I shows the data yield and prediction 
for the initial preselection and for the single- and double- 
tagged samples and the expected Higgs boson signal for 
a Higgs boson of mass 115 GeV/c^. 

The backgrounds from ti, M^-|-jets, Z-|-jets and diboson 
processes are estimated using simulated events. These re- 
sults are validated using the EW enriched control sample. 
The shape of the multijet background is taken from the 
multijet modeling sample, and the normalization is de- 
termined by forcing the number of multijet events plus 
the number of SM predicted background events to equal 
the data yield in the preselection sample. This proce- 
dure is validated by comparing the prediction with the 
multijet enriched sample. 

The signal-to-background separation is then further 
improved using boosted decision trees (BDTs) . For each 
TTiH considered, a multijet boosted decision tree (MJ 
BDT) with 23 input variables is trained on Higgs signal 
and multijet backgrounds. Events which have a MJ BDT 
output greater than 0.6 are retained. These events are 
then input to a second BDT (SM BDT) trained on the 
remaining backgrounds and Higgs signal events using the 
same 23 variables input to the MJ BDT. Fig. [60] shows 
the BDT outputs for the data, predicted background and 
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FIG. 60 The outputs for the (a) MJ BDT, (b) the SM BDT 
for the single-tagged channel and (c) the SM BDT for the 
double-tagged channel for the DO VH — >■ $j.bb analysis. 



signal. The agreement between data and prediction is 
good. Limits are extracted by fitting signal and back- 
ground SM BDT outputs to the data distribution using 
the same modified frequentist algorithm as was used for 
the previously described DO results. Fig. [61] shows the 
Higgs cross section limits from the DO analysis. 

The systematic uncertainties for the DO result are dom- 
inated by similar sources as for the CDF analysis. The 
DO jet energy calibration systematic is < 10%, and the 
background normalization systematic varies between 6% 
- 20%. DO also reports additional systematic uncertain- 
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Sample 


ZH 


WH 


W+jets 


Z+jets 


Top 


VV 


Multijet 


Total background 


Observed 


Preselection 


13.73 ± 1.37 


11.64 ± 1.17 


19 069 


9432 


1216 


1112 


1196 


32 025 ± 4037 


31718 


Single-tagged 


4.16 ± 0.42 


3.60 ± 0.37 


802 


439 


404 


60 


125 


1830 ± 255 


1712 


Double-tagged 


4.66 ± 0.58 


4.00 ± 0.50 


191 


124 


199 


24 


< 8 


538 ± 87 


514 



TABLE XXXVIII The predicted and observed yields for the DO VH — j- ]pj.hb analysis for the untagged control sample and the 
single- and double-tagged analysis samples after applying the multijet BDT selection. 
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FIG. 61 The expected and observed 90% CL upper bounds 
on the Higgs boson production cross section in the DO VH — > 
lfij,hb search. 

ties from luminosity, trigger and identification efficiencies 
and ^-tagging efficiency. These sources have an uncer- 
tainty of roughly 5% for each. 

4. r+jets and rr+jets Final States 

The searches described earlier for the processes WH — 
ivbb and ZH — >• Ubb required £ = e, fi, but correspond- 
ing channels exist with £ — t leptons in the final state. In 
addition, B{H tt) « 7% in th e SM, the most impor- 
tant mode after H — > bb. D0 |l5Cll | has published searches 
for Higgs production in the modes VH — rbb t^v 
and VH rrjj -> ^ivVrThjj, pp VVjj Hjj -> 
TTjj fivi^rThjj, and gg -> Hjj rrjj HWrThjj 
using a data sample corresponding to 1.0 fb^^. Here 
represents a r decay to hadrons. 

The primary signal contribution to the Tvvbb final 
state comes from WH — >■ ri/bb, but roughly 7% of the 
final signal acceptance comes from ZH — > rrbb in which 
one r is not identified. The identified r is required to 
decay to hadrons 58]. Candidate events are initially se- 
lected using a trigger requiring a high px jet and large 
$]rp. The offline selection requires a hadronic r decay, 
with a minimum transverse energy from 12 GeV/c to 20 
GeV/c, with the range depending in which of two decay 
topologies the r lepton is reconstructed. The topolo- 
gies correspond roughly to single charged- hadron r lep- 
ton decays (type 1) and single charged- hadron with mul- 
tiple neutral-hadron decays (type 2). In addition to the 



T lepton candidate, the basic offline selection requires 
]^rp > 30 GeV, the presence of at least two candidate 
6-jets with pr > 20 GeV/c and \r]\ < 2.5, < 3 jets 
with pt GeV/c, and a veto of events with a high-p^ 
electron or muon to keep orthogonality with the other 
searches. Events are also required to pass /S.(j)($rp ^"fj,) < 
7r/2, Ht < 200 GeV, 50 < Mjj < 200 GeV/c^ and 
A(j){T,$T) < 0.02(7r - 2)(^j. - 30) -I- 2. In addition, for 
events with type 2 tau decay, the transverse mass of the 
f and ^ must satisfy Mr < 80 GeV/c^. 

The TTjj final state search has signal contributions 
from a variety of processes: ZH — > rrbb, ZH — >■ qqbb (in 
which H — T' TT, VBF pp — >■ VVjj — > Hjj — > rrjj — > 
livUrThjj, and gluon fusion gg Hjj -> tt^j ^■ 
fivVrThjj- Unlike other low mass channels, some of these 
processes have no 6-quarks in the decay chain. Conse- 
quently, no 6-jet identification is required for this search. 
One of the t leptons is required to decay to hadrons, and 
the other t lepton is required to decay to a muon. Here 
three topologies are used. Two of these are the same 
as in the Tubb search, and the third corresponds roughly 
to T lepton decays to three charged hadrons. Candidate 
events are initially selected by requiring at least one trig- 
ger from a set of single muon and muon plus jet trig- 
gers. The initial offline selection requires a muon with 
Pt > 12 GeV/c and rj < 2.0, a hadronic t lepton de- 
cay candidate with a minimum transverse energy from 
15 GeV to 20 GeV (again depending on topology) and 
at least two jets with pT > 15 GeV/c and ri\ < 2.5. No 
6-jet identification is used. Signal and background sepa- 
ration is then improved using a set of neutral networks. 
A separate neutral network is trained for each of the four 
signal sources {WH, ZH [H bb), ZH [H tt) 
and VBF), paired with each of the main backgrounds 
(W^+jets, Z-|-jets, it and multijet). In addition sepa- 
rate networks are trained for low mass (105, 115 and 125 
GeV/c ) and high mass (135, 145 GeV/c ) giving a to- 
tal of 32 networks. Each network uses six or seven input 
variables chosen such that each one individually improves 
the S/ \fB when it is added as the last variable. The same 
choice of input variables is made at all masses for each 
signal and background pairing. 

All backgrounds except that from multijet events in 
which jets are misidentified as t leptons are estimated 
using simulated events processed through a detailed de- 
tector simulation with corrections based on comparison 
of data and simulated control samples. The absence of 
electrons and muons from direct W or Z decay results in 
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Search Channel 


Source 


Ti/bb 


TTjj 


W+lf 


0.5 ±0.0 


5.1 ±0.3 


W+M 


10.9 ±0.3 


0.9 ±0.1 


Z+lf 


< 0.2 


43.8 ±0.6 


Z+M 


0.4 ±0.0 


10.1 ±0.7 


ti 


9.5 ±0.1 


2.8 ±0.0 


Diboson 


U. / It u.u 


9 1-1-00 


Multijet 


1.3 ±0.1 


6.5 ±2.8 


Total 


23.3 ±0.4 


71.2 ±3.0 


Data 


13 


58 


Signal 


0.216 


0.293 



TABLE XXXIX Expected and observed yields for the DO 
VH — )■ Tvbb and VH — >■ rrjj searches. Only statistical un- 
certainties are shown. 



a significantly higher multijet background in this search 
than in most. The multijet background is estimated from 
using a control data sample enriched in r-like hadronic 
jet events. Systematic uncertainty sources include trigger 
efficiency, lepton and jet identification efficiencies, jet en- 
ergy calibration and production cross section uncertain- 
ties for backgrounds determined using simulated events. 

As for the other analyses, systematic uncertainties are 
divided into those which affect only normalization and 
tliose wliich affect the shape. Sources affecting the nor- 
malization include: the integrated luminosity, the trigger 
efficiency, muon identification efficiency, r lepton identifi- 
cation efficiency, jet identification efficiency, the r lepton 
energy calibration and background cross sections. The 
systematic uncertainties which affect the shape for the 
Tvbb analysis are the jet energy resolutions, the jet en- 
ergy calibration and the 6-tagging efficiencies. For the 
TTjj analysis only the multijet background had a shape 
dependence. The systematic uncertainties, expressed as 
a fraction of the related source, range between 3% and 
30%. For example, the W^±hf cross section systematic 
is 30% of the VF+hf background, but because the back- 
ground is roughly 50% of the total background in the 
Ti/bb channel, the effective systematic is roughly 15%. 

Predicted background and signal levels and observed 
yields for the two searches are shown in Tab. IXXXIXI 
assuming a signal mass of Mh = 115 GeV/c^. The dis- 
tributions or Mjj and the output from the Z±jets NN 
NNz are used as inputs to the limit setting program for 
the Ti^bb and rrjj searches respectively. The distribu- 
tions are shown in Fig. 1621 and the resulting limits are 
shown in Tab. IXLI 



5. H^WW* 

For mH ^ 135 GeV/c^, the dominant Higgs decay is 
to W boson pairs. Below the W pair production thresh- 
old, one of the W bosons will be off-shell which leads to 
lower pt leptons from the W* decay. The most sensitive 
high mass Higgs channel is gg ^ H ^ WW'^*^ — iviv, 




Diboson 
— Higgs x 30 



150 200 250 300 

Dijet mass (GeV) 




0.6 0.8 1 

NN output 



FIG. 62 The dijet invariant mass distribution for the DO 
VH — >■ Tvbb search, and the NN output for the VH — s- rrjj 
search. 



where the two charged leptons in the final state of are of 
opposite charge. In addition to the dominant ggH pro- 
cess, other Higgs production processes that contribute at 
the Tevatron are VBF and VH where V = {W, Z). 

Earl y Tevatron Run II searches for a high mass Higgs 
boson jl5l| . |l52| focused mainly on the ggH process 
and exploiting the angular correlation between final state 
charged leptons due to the scalar (spi n-0) nature of the 
SM Higgs. In a more recent analysis [l53j |. CDF used a 
NN technique to search for a Higgs boson signal in dilep- 
ton + $!rp events having either zero or one reconstructed 
jet. Several kinematic variables, including the results of 
matrix element calculations that combined charged lep- 
ton and information, were used as inputs to the NN. 

The CDF and DO searches for H WW'-*^ in [fsi ] 
and [l55j |. respectively, significantly increase their sensi- 
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TABLE XL Expected and observed 95% CL signal cross- 
section upper limits for the DO VH rvbh and VH — >■ Trjj 
searches. The limits are given as a ratio of the cross-section 
limit to the SM cross-section prediction. 



tivities compared to previous results through the use of 
additional data, topologies arising from VBF and VH 
channels, and analysis improvements such as increased 
charged lepton acceptance (DO). 

In the DO analysis [155] , an integrated luminosity of 5.4 
fb^^ is used to search for H — lyVF'*' in events with two 
oppositely-charged leptons e"'"e~, e^fj,^, or Elec- 
trons are required to have \r]\ < 2.5 (< 2.0 in the e+e^ 
channel) and > 15 GeV. Muons are required to have 
\r]\ < 2.0 and pt^ > 10 GeV/c (in the fi^fi^ channel, one 
of the two muons is required to have pt^ > 20 GeV/c). 
In addition, the dilepton invariant mass is required to 
exceed 15 GeV/c^. Reconstructed jets are required to 
have E^r^ > 15 GeV and jryj < 2.4, however no jet-based 
event selection is applied since the number of jets in the 
event is used as input to an NN to help discriminate sig- 
nal from background. The dilepton invariant mass after 
this "preselection" is shown in Fig. [Mia). 

Additional cuts are made to suppress Z/j* production. 
These include requiring > 20 GeV (> 25 GeV in the 
fj-^ fj-^ channel), high quality measurement, a mini- 
mum W transverse mass to be > 20 GeV (> 30 GeV in 
the e+e~ channel), and azimuthal opening angle between 
the two leptons A(j){i,£) < 2.0 rad. 

To improve the separation between signal and back- 
ground, a NN optimized for several mn values between 
115 GeV/c^ and 200 GeV/c^ is used in each of the three 
channels. Several discriminant variables are used as in- 
puts to the NN: the transverse momenta of the leptons, 
a variable indicating the quality of the leptons' identifi- 
cation, the transverse momentum and invariant mass of 
the dilepton system, minimum W transverse mass, $!rp, 
$j, quality, A0(^,^), A0(£i,^j,), A(f>{£2, $Jt): t^ie num- 
ber of identified jets, and the scalar sum of the transverse 
momenta of the jets. 

Fig. 123] shows the agreement between data and Monte 
Carlo simulation after final selection for (b) the A(j){£, £) 
angle and (c) the neural network output. The expected 
and observed event yields are shown in Tab. IXLII 

No significant excess of signal-like events is observed 
for any test value of itih after the final selection. The 
NN output distributions are used to set upper limits on 
the SM Higgs boson production cross section. Fig. IMl a) 
shows a comparison of the NN distribution between 
background-subtracted data and the expected signal for 



rriH = 165 GeV/c^ hypothesis. Fig. [Mlb) shows the 
expected and observed upper limits as a ratio to the ex- 
pected SM cross section. Assuming niH = 165 GeV/c^, 
the observed (expected) upper limit at 95% CL on Higgs 
boson production is a factor of 1.55 (1.36) times the SM 
cross section. 

In the CDF analysis [154| , an integrated luminosity of 
4.8 fb~"^ is used to search for H — >• VFVF^*^ in events with 
either an opposite sign (OS) or same-sign (SS) charged 
lepton pair. This is an inclusive search that expands the 
signal acceptance by 50% for mn = 160 GeV/c^ com- 
pared to searching for only the g pH production process 
as published previously by CDF [l53j . 

The CDF analysis uses physics objects identified as 
jets, electrons, and muons as well as the $Jrp in events. 
The search is based on the requirement that events con- 
tain two charged leptons resulting from the decays of the 
final-state vector bosons which are OS except in the case 
of the VH channel where they can be SS. At least one 
charged lepton is required to match the lepton found in 
the trigger and have Et{pt) > 20 GeV (GeV/c) for elec- 
trons (muons) . The second charged lepton is required to 
have Et{pt) > 10 GeV (GeV/c) except in events with 
same charge leptons, where both leptons are required to 
have Et{pt) > 20 GeV (GeV/c). Requirements on the 
event $Jrp indicative of the presence of neutrinos from W 
boson decay are made in opposite-charge dilepton events. 
Backgrounds due to Drell-Yan and heavy flavor are sup- 
pressed by requiring that the invariant mass of the lepton 
pair be greater than 16 GeV/c^. The Higgs boson signa- 
ture can also involve jets of hadrons produced from the 
decay of one of the vector bosons in the VH process, for- 
ward quarks in the VBF process, or from the radiation 
of gluons. 

In order to increase the sensitivity to the various Higgs 
production processes in the SM, candidate events are sub- 
divided into six analysis channels based on jet multiplic- 
ity, lepton categories, and lepton charge combinations. 
Five of the channels have signatures with OS leptons and 
the other is for SS leptons. 

Discrimination of signal from background is based on 
NNs trained in each analysis channel and at each of 
14 hypothesized niH values in the range 110 < mn < 
200 GeV/c. The NN inputs are based on kinematic 
quantities selected to exploit features such as the spin 
correlation between the W bosons in Higgs boson decay, 
the presence of large $!rp from the neutrinos; the trans- 
verse mass of the Higgs boson (calculated using the lep- 
tons' four- momenta and $!rp vector); and the modest total 
energy of the Higgs boson decay products compared to it 
decay. In the zero jet categories, we additionally classify 
events by evaluating the observed kinematic configura- 
tion in a likelihood ratio of the signal probability density 
divided by the sum of the signal and background prob- 
ability densities. These probability densities are deter- 
mined from LO matrix elem ent calculations of the cross 
sections of each process jl53l | (see also Section lTlI.E.4l for 
the case of SM WW). 
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FIG. 63 For the DO H — > W^W^'*' analysis [155| ]. the (a) dilepton invariant mass after preselection; (b) /S.(f){l,t) angle after 
final selection; and (c) neural network output after final selection. The signal is shown for mH=165 GeV. The systematic 
uncertainty is shown after the fit. 
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TABLE XLI Expecte d an d observed event yields in each channel after preselection and at the final selection for the DO 
H — !> VK'py'^*^ analysis [155| ]. The systematic uncertainty after fitting is shown for all samples at final selection. 



An example NN discriminant distribution for the com- 
bination of all categories is shown in Fig. [SSJ where signal 
and background expectations for a 160 GeV/c^ Higgs bo- 
son are compared to the observed data. 

No significant excess of events beyond SM background 
expectations in the NN discriminant are observed. The 
95% confidence limits on an, expressed as a ratio to the 
expected SM rate as a function of ttih-, are determined 
from the data. The median expected and observed upper 
limits on au as a function of rriH is shown in Fig. 1661 

The expected and observed upper limits at 95% CL for 
Higgs boson production cross section expressed as a ratio 
to the cross section predicted by the SM as a function of 
niH for t he C DF and DO searches for H — > M^M^^*) in 
|154| and |l55j |. respectively, is shown in Tab. IXLlIl 



C. Combined Limits 

Results for individual search final states are combined 
by the Tevatron New-Phenomena and Higgs Working 



Group|156!|. The combinations are updated roughly ev- 
ery six months to one year using the most up to date 
preliminary and published results. The combination is 
performed using a joint probability density incorporat- 
ing each channel as an individual result with systematic 
uncertainties and their correlations included. Two algo- 
rithms ar e us ed for the combination. One is a Bayesian 
approach 154 , an d the other is a modified frequentist 
approach 157 -[ll^ . In both algorithms, the impact of 



systematic uncertainties is reduced by treating these as 
nuisance parameters. Results for the 95% CL upper 
bound on the Higgs production cross section from the 
two algorithms agree to within 1 0%. 

The most recent combination [iS^ incorporates pub- 
lished and preliminary results based on samples rang- 
ing from JC dt = 2.1 fb'^ to J C dt = 5.4 ib'^, and 
combined H — >• WW limits are at jl60l |. The cross sec- 
tion upper bounds are a factor of 2.7 higher than the 
SM prediction for rriH — 115 GeV/c^ and production 
of a SM Higgs is excluded at the 95% CL in the region 
163 < fUH < 166 GeV/c . The results are shown in 
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TABLE XLII Expected and observed upper limits at 95% CL for Higgs boson production cross section expressed as a ratio to 
the cross section predicted by the SM for a range of test Higgs boson masses. 
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FIG. 64 For the DO H ^ WW^*> analysis ^5^, the (a) 
data after subtracting the fitted background (points) and SM 
signal expectation (filled histogram) as a function of the NN 
output for niH = 165 GeV /c^ . Also shown is the ±1 standard 
deviation (s.d.) band on the total background after fitting, 
(b) Upper limit on Higgs boson production cross section at 
95% CL expressed as a ratio to the SM cross section. The 
one and two s.d. bands around the curve corresponding to 
the background-only hypothesis are also shown. 
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FIG. 65 The combined distribution of NN scores for 
grounds and a Mh = 160 GeV/c^ Higgs boson compared to 
the observed data shown with statistical uncertainties for the 
WW"'*) 



CDF H WW^*' analysis [15J]. The Higgs boson distribu- 
tion is normalized to ten times the SM expectation. 
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FIG. 66 Expected and observed upper limits at the 95% CL 
on an presented as a ratio to the predicted SM values as a 
function of Mh for the CDF H -> WW**' analysis [HI. The 
dashed line represents the median expected limits, the green 
and yellow bands the estimated one and two sigma probability 
bands for the distribution of expectations, and the solid line 
the observed limit. 
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FIG. 67 The 95% CL upper bounds on the Higgs production 
cross section divided by the predicted cross section. This 
result includes both published and preliminary results from 
CDF and DO. 



Fig. [S71 expressed as a ratio of the 95% CL upper bound 
cross section divided by the SM prediction. 



VI. SUMMARY AND OUTLOOK 

The SM has been tested extensively over the last forty 
years. The only significant deviation from its predictions 
is the existence of neutrino mass. Ongoing studies of the 
SM are being carried out by the Tevatron experiments, 
CDF and DO. This paper presents results of these tests 
and discusses prospects for the remainder of the Tevatron 
running and the transition from the Tevatron to the LHC. 

CDF and DO have presented a wide variety of results 
relating to the EW sector. Among the results are ob- 
servation of all SM diboson processes involving W, Z 
and 7, including for the first time WZ production, ob- 
servation of EW production of single top quarks and the 
corresponding measurement of \Vtb\, world's best mea- 
surements of the top quark and W boson masses allowing 
improved indirect constraints on the Higgs boson mass, 
and ongoing searches for the Higgs boson which have be- 
gun to exclude mass values outside of the existing limits. 
In addition to testing the SM, some of these are bench- 
mark measurements eventually to be compared with cor- 
responding results from the LHC. Others, like the top 
quark and W boson masses, are likely to be legacy mea- 
surements from the Tevatron. 

With additional luminosity still to come for CDF and 
DO, many of these measurements will significantly im- 
prove before the end of Run II at the Tevatron. This is 
an exciting time in particle physics as the field transitions 
from the era of the Tevatron to an era of the LHC as the 
machine to probe physics at the energy frontier. 



The many beautiful experimental measurements sum- 
marized in this paper are the result of the work of thou- 
sands of people from the Fermilab accelerator division, 
CDF and DO. The authors wish to thank those involved 
in making these possible. In addition, the authors thank 
Andrew Askew, Tom Diehl, Reinhard Schwienhorst and 
Junjie Zhu for useful discussion and clarification. 
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